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ABSTRACT 


VLF wave-wave interaction experiments have been carried out by in- 
jecting various forms of VLF pulses into the magnetosphere from a 21.2 
Km dipole antenna at Siple, Antarctica. The data are collected at the 
Siple conjugate, near Roberval, Quebec. The injected signals propagate 
along a geomagnetic field line and often interact strongly with ener- 
getic electrons trapped in the radiation belts near the equator. Sig- 
nals may be amplified and trigger emissions. These signals may then 
interact with one another through these energetic electrons. Suppres- 
sion of signal growth and energy coupling among the signals are commonly 
observed . 

This report is divided into three parts. In the first part, simu- 
lations of VLF pulses propagating in the magnetosphere are carried out. 
The results show that the leading edge of a VLF pulse is distorted suf- 
ficiently when it arrives at the interaction region to alter the details 
of the wave-particle interaction process in the magnetosphere. It is 
also found that the calculated group delay at the carrier frequency agrees 
closely with the travel time of the center of a VLF pulse. The center is 
defined as the center of the interval within which the signal exceeds 30 
percent of the peak value. In addition, a digital equalizer is designed 
to compensate for the distortions developed in a VLF pulse that propagates 
through the magnetosphere. Computer simulations show that this device 
can completely compensate for the distortions. 

In the second part, it is found for the first time that a 10 ms gap 
in a triggering wave can induce emission, which may then interact 
with the post-gap signals. Suppression and entrainment between the 
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gap-induced emissions and the post-gap signals are often observed. The 
Phenomenon of gap-induced emission is explained by the hypothesis that 
an emission can develop within 10 msec at the end of a triggering wave. 
Experimental data show that, as the gap is opened up in time, we can 

"scan" the development of emissions. One of the results confirms an 
earlier observation that emissions always start with a small rise in 
frequency. It is also found that a 70 ms gap in a trigger wave is long 
enough to allow falling emission fully developed. During these studies, 
a new feature of wave-wave interaction was observed. The post-gap sig- 
nals can "capture" a falling emission that is 90 Hz below the transmitter 
signals, causing the emission to overshoot toward the positive frequency 
side. The frequency difference between these two signals appears to be 

an important factor in determining whether or not the interaction can 
occur. 

In the third part, sideband triggering is reported for the first 
time. Sidebands are generated by frequency-shift-key (FSI0 modulation 
techniques. Sidebands generated by short modulation periods can trigger 
emissions, indicating that the average electron "memory" time in the mag- 
netosphere is at least as long as twice the modulation period. The cor- 
responding length of the electron interaction region is estimated to lie 
between 2000 and 4000 km. Furthermore, it is found that two sidebands 
with a frequency difference of 50 Hz or less mutually interact. The in- 
teraction is attributed to the overlap of the two Vj| ranges in which 
the electrons can be organized by the corresponding sidebands. It is 
believed that the V ( range is related to the wave intensity in the in- 
teraction region. Using a homogeneous model, we have found the wave in- 
tensity in the interaction region to be 2.5 ~ 10 my, in reasonable 
agreement with satellite measurements. 



iv 


ACKNOWLEDGMENTS 


I would like to thank Professor R. A. Helliwell for his guidance of 
this research. Professor Helliwell has impressed me with his professional 
integrity and his constructive criticism of new ideas . His insight and 
comments have improved both the scientific content and the literary style 

of this report . 

I owe a great deal to Dr. T. F. Bell for revising the manuscript and 
his suggestions and comments at various stages of the work. I am grate- 
ful to Professors D. L. Carpenter and J. K. Katsufrakis and Dr. C.G. Park 
for many suggestions and comments . Most of the data presented here are 
related to the transmitter at Siple Station, Antarctica. My sincere 
thanks go to the various people who went to Siple Station. I am grate 
ful to Mr. J. W. Yarbrough for making all the spectrograms in this report. 
Informal discussions with my colleagues in the Radioscience laboratory, 
Stanford University, especially Mr. E. C. Ha with whom I have shared the 
same office for three years, have contributed significantly to this work. 

I would like to thank Professor M. Chodorow for his careful reading 
of this work. I want to thank all persons involved in the final prepara- 
tion of this manuscript, in particular , ingaUf, for her excellent typing. 

I am grateful to my parents for their moral support and to my wife, 
Bi-Hua, for her patience and encouragement throughout the course of this 

work . 

This research was supported in part by the Division of Polar Programs 
of the National Science Foundation, in part by the Atmospheric Science 
Section of the National Science Foundation, and in part by the National 
Aeronautics and Space Administration. 


v 


CONTENTS 


Page 

I . INTRODUCTION x 

A. Magnetosphere 

B. Whistler-Mode Signals and Cyclotron Resonance 4 

C. Background of the Experiments 5 

D. Outline of the Thesis 10 

E. Contributions of the Present Work 

1. Pulse Propagation in the Magnetosphere 12 

2. Gap-Induced Emissions 13 

3. Sideband Triggerings . 14 

II. VLF PULSE PROPAGATION IN THE MAGNETOSPHERE 

A. Introduction 

B. Theory 

1. Homogeneous and Nondispersive Medium 21 

2. Dispersive Medium . 22 

C. Propagation in a Homogeneous Model 27 

D. Pulse Propagation in the Magnetosphere o 34 

E. Discussion 3g 

F. Conclusion ........ 48 

III. GAP INDUCED EMISSIONS 51 

A. Introduction 

B. Examples of Gap-Induced Emissions 53 

C. An Interpretation of Gap— Induced Emissions ....... 64 

D. GATR Transmitter Program 

E. Summary and Recommendations . 7g 

IV. SIDEBAND TRIGGERING S1 

% 


me. 


vii 



CONTENTS (Cont) 


Page 

A. Introduction 81 

B. A Transmitter Program for Generation of Sidebands .... 82 

C. Estimation of the Length of the Interaction 

Region in the Magnetosphere 9° 

D. Sideband Triggering and Features of Coupling 101 

E. Discussion -*-08 

F. LISP Program 3-20 

G. Summary and Conclusion 131 

V. SUMMARY AND RECOMMENDATIONS 133 

A. Introduction *- 33 

B . Summary of the Present Work I 33 

1. Pulse Propagation in the Magnetosphere . 133 

2. Gap-Induced Emissions 124 

3. Sideband Triggering 133 

C. Suggestions for Future Work 137 

1. pulse Propagation in the Magnetosphere 137 

2. FNWWI 137 

3. SMI . . 137 

Appendix A. A COMPUTER CODE SIMULATING VLF PROPAGATION 

IN THE MAGNETOSPHERE 139 

Appendix B. VLF WAVE-PARTICLE INTERACTION IN THE 

MAGNETOSPHERE ........ 145 

Appendix C. DERIVATION OF SIDEBAND REPRESENTATION OF 

FSK WAVES 173 

Appendix D. MISCELLANEOUS EXAMPLES 183 

BIBLIOGRAPHY 199 


viii 


ILLUSTRATIONS 


Figure Page 


1.1 

A sketch of the noon-midnight meridian section 
of a model magnetosphere showing the plasmasphere 
(shaded area inside the plasmapause) in relation 
to other magnetosphere features 

... 2 

1.2 

A sketch showing propagation of transmitter sig- 
nals between Siple station ana Roberval 

... 7 

1.3 

Typical examples of TO wave amplifications and 
triggering emissions 

... 8 

2.1 

Flow chart of the method of investigating VLF 

pulse propagation in the magnetosphere ... 

... 25 

2.2 

Demonstration of the sampling theorem . 

... 26 

2.3 

Decomposition of a circularly polarized signal 

into two linearly polarized signals 

. . . 28 

2.4 

A plot of group velocity vs frequency of a TO 

signal propagating in a homogeneous magnetoplasma 

in which f =180 kHz and f„ - 13 kHz 

N H 

... 30 

2.5 

Propagation of a 50 ms pulse at f = 3.25 kHz 
through the homogeneous magnetoplasma 

... 31 

2.6 

Propagation of six 20 ms pulses at 1.75, 2.25, 

2.75, 3.25, 3.75, and 4.25 kHz through the homo- 
geneous magnetoplasma for 25,000 km ......... 

... 33 

2.7 

Propagation of 10 ms pulses at 5.065 kHz through 
a duct at L = 4 in a DE-1 model magnetosphere . . . 

... 37 

2.8 

Propagation of 10 ms pulses at 7.065 kHz in the 
model magnetosphere 

. . . 39 

2.9 

Propagation of 10 ms pulses at 6.065 kHz in the 
model magnetosphere 

... 40 

2.10 

Propagation of 10 ms pulses at 4.065 kHz in the 
model magnetosphere 

• 0 9 41 

2.11 

Propagation of 10 ms pulses at 3 .065 kHz in the 
model magnetosphere 

... 42 


2.12 Propagation of a 100 ms pulse at 6.6 kHz with a 

phase discontinuity through a duct at L = 3.1 47 


ix 


ILLUSTRATIONS (Cont) 


Figure 

Page 


2.13 

Propagation of a pre-processed 100 ms pulse (with a 
phase reversal) at 6 .6 kHz through a duct at L = 3.1 

• • • • 49 

3.1 

Demonstration of the idea of how a phase ■ reversing 
signal affects the phase-bunched current . 

• » ♦ # 54 

3.2 

Demonstration of how a phase reversal is achieved 

• * • * 56 

3.3 

Two examples of typical gap-induced emissions . . . 

. . . . 57 

3.4 

Examples of gap-induced emissions . . . 

.... 60 

3.5 

Examples showing wave -wave interactions between 
gap-induced emissions and the post-gap signals . . . 

• * • • 61 

3.6 

Electron motion in an 0 wave . . . 

. . . . 66 

3.7 

Electron motion in a rt wave .... 

• . . . 67 

3.8 

Deformation of the trapping domain during the gaps . 

• • • • 69 

3.9 

Deformation of the trapping domain during a 10 
msec gap in an inhomogeneous model . . . 

. . . . 71 

3.10 

The features of gap-induced emissions vs various 
gap sizes 

. . . . 77 

4.1 

Waveform and f-t display of an FSK signal 

. . . . 84 

4.2 

Part of the transmitter program described in 
Section B 

. . . . 85 

4.3 

Two different f-t displays of an FSK signal 

• • • ss 

4.4 

Histogram of the number of the FSK waves that 
trigger emissions vs the values of T 

s 

... 93 

4.5 

Examples of the FSK waves from which Fig. 4.4 is 
constructed 

... 97 

4.6 

Illustration of the apparent interaction time (AIT) 

... 99 

4.7 

A typical example of sideband triggerings . . . 

. . . 102 

4.8 

More examples of sideband triggerings 

... 104 

4.9 

An example of sideband couplings .... 

... 107 

4.10 

An electron motion in two waves . . . 

... Ill 


x 



ILLUSTRATIONS (Cont) 


Figure Page 


4.11 Demonstration of the ranges of V |( over which 
the electrons are organized by two waves with 

a large frequency difference 114 

4.12 Plots of f(Q!) = sin CL cos^ V a as a function 

of d for various values of m and p 117 

4.13 Plots of r(m,p) for the case where CL = 7° 

and Cty =87° 119 

4.14 A sketch of the LISP program that consists of 

FSK waves with t = 10 ms 121 

s 

4.15 The spectrum of an FSK wave with Af = 100 Hz 
and t s = 10 ms and that with Af = 200 Hz 

and t s - 10 ms 122 

4.16 An example of the LISP program 123 

4.17 Expanded display of the dynamic spectrum in 

the eighth minute in Fig. 4.16 125 

4.18 Expanded time scale of the transition at the 

second minute in Fig. 4.16 .. 127 

4.19 Expanded time scale of the transition at the 

fourth minute in Fig. 4.16 128 

4.20 Expanded time scale of the transition at the 

eighth minute in Fig. 4.16 129 

4.21 Expanded time scale of the transition at the 

tenth minute in Fig. 4,16 190 

B.l Configuration of E and B fields of a 

whistler mode wave and velocity components 

V l( and V x of an electron in the lab frame 147 

B .2 The phase-plane loci of the vertices of the 

parabolas defined by Eq. (B.4b) 149 

B.3 Three surfaces in the phase plane specified 
respectively by three constants X; X > 0, 

X = 0, and X<0 150 

B.4 Trapping contours in the 3-D phase space 151 

B.5 The angle vs wave amplitude, assuming 

electron pitch angle to be about 35° to 38° 153 

xi 


ILLUSTRATIONS (Cont) 


Figure 


B.6 The projection of the 3-D trapping contours shown 
in Fig. B.4 onto a (V|| — \Jr ) plane ....... 

B.7 The total trapping width AV tu + AV +T and the 
value of AV^ - AV tu vs wave intensity .... 

B.S The calculated trapping width AV t by Eq. (B.10) 
and the value of AV tu + AV tL vs wave intensity , 

B.9 An analog of an electron motion in a TO wave to a 
pendulum motion in the gravitational field . 

B.10 Plot of the value of K(sin ijr /2) in Eq. (B.14b) 
vs * . . . R 


B.ll Illustration of phase bunching mechanism . . 

B.12 Configuration of j j of Eqs . (B.17a> and 
(B .17b) ..... T 


B * 13 Variation of a trapped region at the end of 

triggering wave in a homogeneous magnetoplasma 

C.l Plot of the functions M^t), M rt (t) of Eq. 
(C.2b) vs time 


C .2 


The corresponding spectrum 

m (f) e i2Tt: fTc 
o 


m. 


r s in the frequency domain 


(f) of M 0 ( t ) and 


C .3 


The corresponding spectrum m_.(f) of M.(t) 
m^(f) e 1 T s i n the frequency domain . . . 


and 


C.4 Plots of A + (f) of Eq. (C.5c) for two O-wave 
cases: t s = 10 ms, Af = 100 Hz and T = 20 
ms, Af = 100 Hz ..... s 


C ' 5 Plots °f A + (f) of Eq. (C.5c) for two x-wave 

cases: t s = 10 ms , Af = 150 Hz and r. = 20 
ms, Af = 125 Hz 

D.l Display of a ’’rake-spectrum" wave at a high time 

lesolution and that at a high frequency resolution 

D.2 An example of multiple triggering; a typical 

feature of triggering at equinox ......... 

D.3 Demonstration of a part of the C0BA program 
transmitted around equinox ...... 


Page 

. . 155 

. . 156 

. . 158 

. . 160 

, . 163 

. 163 

. 170 

. 170 

. 174 

. 177 

. 178 

. 179 

. 180 
. 184 

. 186 
. 186 


xii 


ILLUSTRATIONS (Cont) 


Figure Page 


D.4 An example of a "rake-spectrum” wave 188 

D . 5 Another example of a "rake-spectrum" wave 189 

D .6 An example of the variations of the saturation 

amplitude of FSK waves as Af increases 192 

D.7 The plot of the amplitude vs Af of the FSK 

waves illustrated in Fig. D.6 193 

D.8 The plot of the amplitude vs Af of the same 
data as that in Fig. D.7 but from a filter 

with a wider bandwidth, 680 Hz 195 

D.9 The plot of averaged amplitude vs Af of the 
FSK waves over four sets of data samples taken 

within an hour . 196 

TABLES 

Number 

3.1 Averaged initial growth rate, post-gap growth 
rate, and saturation level over 52 examples 
taken on April 28 and May 13, 1975, for both 

0-waves and 63 

4.1 The relative phase between successive pulses 
at a particular carrier frequency for various 

FSK waves * • 86 


xiii 



Chapter I 


INTRODUCTION 

This study is concerned with VLF wave-wave interaction (WWI ) in the 
magnetosphere. In this chapter, we give a brief introduction on this 
topic and outline the contributions of the present work. 


A . Magnetosphere 

The magnetosphere is created by the interaction between the solar 
wind and the geomagnetic fields. The solar wind is a plasma consisting 
mainly of electrons and protons. The solar wind flow is supersonic. At 
the earth, the solar wind encounters the earth's magnetic field. The 
result is the bow-shock wave sunward of the planet, where the solar wind 
is decelerated and deflected. The solar wind drags the geomagnetic field 
toward the midnight direction. The magnetotail is formed. Most of the 
solar wind flows past the earth's protective magnetic field, but some 
leaks through, filling parts of the magnetic cavity with hot plasma. 
This magnetic cavity, distorted by solar wind flow and filled with ten- 
uous plasma of varying temperatures, is the magnetosphere [e.g., Chapman, 
1977]. A sketch of the magnetosphere is shown in Fig. 1.1. 

The region of our interest is the so-called "plasmasphere that, 
roughly speaking, corotates with the earth. The plasmasphere is popu- 
lated by cold plasma, mainly electrons and protons in the energy range 
of a few electron volts or less. The equatorial density could be as 
high as 10 2 to 10 3 electrons per cm 3 . This equatorial density drops by 
one to two orders of magnitude within 0.5 earth radii (R^) or less at 
the boundary between the outer magnetosphere and the plasmasphere 
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BOW SHOCK 



Fig. 1.1. A SKETCH OF THE NOON-MIDNIGHT MERIDIAN SECTION OF A 
MODEL MAGNETOSPHERE SHOWING THE PLASMASPHERE (SHADED AREA INSIDE 
THE PLASMAPAUSE) IN RELATION TO OTHER MAGNETOSPHERE FEATURES . 
(Taken from Ho [1974]) 


[Carpenter, 1962, 1966]. The boundary, discovered by Carpenter [1962, 
1963], is called plasmapause. The location of the plasmapause varies 
from 2R^ after the onset of a major substorm to as much as 7R^ after 
a long quiet period. The plasmapause is generally near ~4R under 
moderately disturbed conditions [e.g., Carpenter, 1966]. 

The cold plasma in the plasmasphere is created mainly by the sun. 
Photons from the sun produce ionization in the F-region of the ionosphere 
in the dayside sector. This cold plasma from the ionosphere can diffuse 
upward along the geomagnetic field lines into the plasmasphere surround- 
ing the earth. At night, the cold plasma "stored" in the plasmasphere 
flows downward into the ionosphere to maintain the night time F-region 
[Park, 1970]. 

Superimposed upon the cold plasma is a tenuous high energy component 
with energy in the keV or higher range. Most of these high energy par- 
ticles are trapped in the plasmasphere for reasonably long periods . The 

electron density of the energetic component can be as low as a few tenths 

3 

of an electron per cm . An electron moving among the magneto spheric 
protons is practically collisionless [e.g., Ratcliffe, 1972] and moves 
in a helical path around a geomagnetic field line and bounces back and 
forth between two mirror points ; one in the northern hemisphere and the 
other in the southern hemisphere. Except for gyrating in the opposite 
sense and taking a longer period than electrons, a proton also forms a 
helical path along a geomagnetic field line and bounces back and forth 
between its mirror points. Superimposed upon the helical motion is a 
slow longitudinal drift around the earth, eastward for electrons and 
westward for protons. The oppositely directed drifts of protons and 
electrons constitute a ring current around the earth, flowing westward 
[e.g., Roederer, 1970], 
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magnetosphere? How 


How does this high energy plasma get into the 
do the particles diffuse within the magnetosphere? What are the leakage 
mechanisms of these particles? All of these questions have not been yet 
completely answered. One possibility may begin with solar wind particles 
leaking across the magnetopause, which separates the magnetosphere from 
the outer space, into a boundary layer called the plasma mantle that 
sheaths the magnetosphere just inside the boundary. These particles are 
accelerated along the open magnetic field lines and enter the plasmasphere 
from the polar cusp. Once arriving at the plasmasphere, these particles 
may move to other field lines through cross-L diffusion [e .g ., Hess , 1968] . 
The cross-L diffusion may come from convective E-field or asymmetrical 
distortions of geomagnetic field by, say, solar storm. The loss mecha- 
nism is believed to be pitch angle diffusion either due to cyclotron 
resonance or due to bounce motion resonance with external excitations 
[Robert , 1969; Tverskoy, 1969]. For further information on these topics, 
the reader is referred to Hess [1968], Roberts [1969], Roederer [1970], 
Akasofu and Chapman [1972]. 


B * Whistl er-Mode Signals and Cyclotron Resonance 

U'his tier-mode (IVH) is one of the modes of propagation of electromag- 
nets waves in a plasma in the presence of a static magnetic field and is 
characterized by the following: 

(1) propagation along magnetic field 

(2) sense of rotation of wave field the same as electrons 

(2) wave frequency less than gyro -frequency and electron 
plasma frequency 


4 


In the magnetosphere, TO waves are, roughly speaking, in the kHz 
range. (We shall use the term "VLF signal" to represent WM signals in 
the magnetosphere hereafter.) in this report, we limit ourselves to the 
ducted TO signals only. Most of the ducted TO waves travel along field- 
aligned enhancements in the magnetosphere. This confines the wave nor- 
mal to small angles from the static magnetic field. Ducts for WM waves 
in the magnetosphere, like optical fiber glass, can guide the signals 
from one hemisphere to the other. The signals are then able to penetrate 
the lower ionosphere and reach the ground . 

One of the interesting properties of TO waves is their low group 
velocity. TO waves also propagate with low phase velocity. A typical 
value of the group velocity near the equator at ~4R in the magneto- 
sphere is about one-tenth of the speed of light . 

Since TO waves and electrons gyrate in the same sense, electrons 
may experience constant wave fields for reasonably long periods when the 
doppler-shif ted frequencies of the waves seen by the electrons equal the 
local gyrofrequency . This is the cyclotron resonance condition. For TO 
cyclotron resonance, the waves and the electrons must travel in opposite 
directions. In the magnetosphere, the electrons that can interact with 
the VLF waves are mainly in the keV range. For constant frequency sig- 
nal, the most effective interaction region is thought to be around the 
equator, where inhomogeneity of the geomagnetic fields is a minimum. 

c ' Background of the Experiments 

Many types of magnet os pheric amplifications and emissions are be- 
lieved to result from cyclotron interaction between TO waves and energetic 
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In the VLF wave injection ex- 


electrons trapped in the radiation belts, 
periraents that we than discuss in this report, VLF pulses are injected 
into the magnetosphere from a 21.2 km dipole ant e„na at Siple, Antarctica. 
The signals propagate along an enhancement duct and reach the equatorial 
region, where they effectively interact with high energy electrons stream- 
ing from the opposite direction. As a result of this interaction, the 
Signals are amplified and trigger emissions which are the "new" signals 
at frequencies other than the parent signals. The amplified signals and 
triggered emissions then travel along the remainder of the duct and ar- 
rive at the conjugate point of Siple near Hoherval, Quebec . The geometry 
of the wave injection experiments is illustrated in Fig. 1.2. 

The spectrogram in pig. !. 3 shows an example of wave amplification 
and emission triggering. The vertical and horizontal axes show the fre- 
quency and time scales, respectively. The blackness indicates the wave 
intensity. A two-second pulse at 3 kHz is injected into the magnetosphere 
fro. Siple. The arrival time of the pulse at Roberval is at t . 0. The 
pulse duration is indicated by the time scale. It is observed that the 
signal is amplified and triggers two rising emissions; one at about 0.5 
sec and the other at the end of the triggering wave. 

in addition to producing wave amplification and triggering emissions, 
these experiments show a variety of interesting and important effects. 

For example, the emissions may dump electrons into the ionosphere Tlnan, 
1977; Dingle, 1978], The precipitated electrons may excite x-ray fe.g., 
Rosenburg et al, 19711, may perturb electron density sufficiently to ’ 
affect sub ionospheric VLF propagation [Helliwell et al, 1973], ,„ d „ay 
modulate the local conductivity sc as to cause the excitation of ULF 
waves [Bell, 1976], It Is important, therefore, to understand the mecha- 
nisms that govern growth of VLF waves and the generation of emissions. 
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number of theories used in explaining wave amplification 

trl “ Ph — ^ Bell and Bune man, 1964; 

Helliwell , 1967; Liemohn, 1967- Dysthe 1971 - d 

, uysthe, 1971; Palmadesso and Schmidt, 

1971; Matsumoto, 1972; Brinca iq 7 o. „ , 

' , Helliwell and Crystal, 1973; Nunn, 

1974; Newman, 1977]. 

in one class of the theories, phase-bunched current is employed 

b ' S " HelliWG11 a " d CrySta1 ’ 1973:i - * —tent frequency signal can 

interact with electrons which are nearlv in » 

nearly m resonance with it. As a 

result of this interaction, the ■■phases" of the electrons are organised 
hy the ware. *et transyerse currents are produced and they radiate a 

flSld ‘ ”* " ^ — « electrons is transferred 

UO ^ and ^ aPPear as bei *S amplified. This process 

of organizing electrons in phases is called the nha n 

ea the phase -bunching process 

0hara ° t6riStia «- tequired for the exchange of energy between the 
" aV ° “ e el8 =‘«>ns is called the bunching time. 

Many VIF wave injection experiments have been designed to study the 
Phase-bunching process in the magnetosphere. Two such experiments have 
eh us to the discovery of two new phenomena: gap-induced emission and 

sideband triggering. Ths gap-induced emissions and post-gap transmitter 
signals often interact. The former can suppress the growth of the Zat- 

ter, and the latter may also suppress the rtmmi 

PP ess the development of the former. This 

-ugory Of wave-wave interaction involving emissions as well as 

amplified injected signais, is one of the w. that will be discussed 

“ thIS reP ”' t - “ ““ ™ — - discussed is the mutual 

interaction among sidebands at various frequencies, sidebands are gen- 

erated by injecting frequency-shift-hey <«> signals into the magneto- 
sphere. Sideband interactions oocur when the frequency apacings between 
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slte and the interaction region is essentia! in the interpretation of 

the results of these wave injection experiments . 

in Chapter III, the phenomenon ol gap-induced emissions will 

a - covered in one such experiment in which a 10 -msec 
presented. It was discovered „ 

nTlri nulse to alter the phase of the sig 
gap was inserted in a one-second pulse 

the urogram ms designed, it was believed that a 10- 
nal At the time the program 

ms gap was too short to affect the wave-particle interaction processes 

and that the min effects should come from the phase alterations . . 

a. j , +v,o Hptails of phase-bunching 
was the first program designed to study the 

mechanisms . 

„„„ alterations are immaterial to the growth 
It turns out that the phas 

iq discovered that a 10-ms gap 
of the post-gap signals. Furthermore, it is 

in a triggering wave can induce emissions. The gap-induced emissions may 
then interact with the post-gap signals. Suppressions and entrainments 
between those signals are often observed. This is one "type" of ~ 
m ve interaction (»I> in the magnetosphere we shall discuss m this 
report. This «!, involving an amplified transmitter signal and an emis 

Sion, will be called a forced-mod./natural-mode wave-wave interaction 

■ is said to be in a natural mode because it is not 

(FNVWI) . An emission is said 

. . , 0nce it has been triggered, it is on its 

driven by external signals . Once 

„ f Qrc ed (or driven) mode 

own. Amplified transmitter signals are ^ ^ 

ic, "mnrced mode and nat-ir»- 
because they are driven by external signals . Forced 

mode" adapted In here are defined by cohsiderlng the 

motion process in the magnetosphere as a whole to be a single system, 
in the interaction region, which is adapted by Helliwell CUTM. 
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In Chapter IV, we shall study a phenomenon called sideband trigger- 
ing. It was discovered unexpectedly during a wave injection experiment 
in which the phase of a one-second signal was alternated regularly to 
investigate the details of the phase-bunching processes in the magneto- 
sphere . Because of the regular modulation of the phase of a signal , 
sidebands are generated. We find that some of the sidebands, just like 
isolated constant frequency signals , can trigger emissions . 

The frequency spacings between the sidebands are controllable by 
changing the modulation periods. It is found that, when the frequency 
spacing is small, suppression or coupling between sidebands is often 
observed . The mutual interactions among sidebands belong to the cate- 
gory of sideband mutual interactions (SMI). 

E . Contributions of the Present Work 

1 . Pulse Propagation in the Magnetosphere 

We have, for the first time, calculated the distortion of a 
VLF pulse (of finite duration) that propagates through the magnetosphere. 
Among the results are the following: 

(1) It is found that the group velocity evaluated at the 
carrier frequency predicts the arrival time of a VLF 
pulse better when the time delay refers to the "center" 
and not to the leading edge of the pulse. 

(2) It is found that the frequency of a distorted pulse is 
not constant over the duration of the pulse . The larg- 
est deviations occur at the front and/or the rear ends 
of the pulse . 
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(3) it is found that the frequency at which a VLF pulse suf- 
fers a minimum distortion at the equator is the "nose" 
frequency which is defined as the frequency at which the 
signal has the minimum time delay [e .g. , Helliwell , 1965] . 

(4) It is found that a discontinuity in the phase of a sinu- 
soidal wave introduced on the ground develops into an 
amplitude gap when the signal arrives' at the equator. 

The gap can be as large as 20 msec for the signal at 3 
kHz propagating along a duct at L = 4. 

(5) A digital equalizer has been designed to compensate the 
phase distortion developed in the magnetosphere. This 
device can be used as a pre- as well as a post-processor. 
Computer simulations show that the distortion can be com- 
pletely compensated by this device. 


2 • Gap-Induced Emissions 

(1) Gap-mduced emissions are reported for the first time. A 
10-ms gap in a VLF wave can induce emissions, suggesting 
that emissions start to develop within a time period 
shorter than 10 msec at the end of the VLF pulse. 

(2) The phenomenon of gap-induced emissions is explained by 
the hypothesis that emissions (natural mode) can begin 
at the end of the triggering signal. 

(3) A transmitter program has been designed to test some of 
the ideas developed from the studies of gap-induced emis- 
sion. One of the results confirms an earlier observation 
that falling emissions always start with a small rise in 
frequency. It is also found that a 70-msec gap is long 

enough to allow falling emissions to become fully devel- 
oped . 
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(4) A new feature of the WWI has been discovered. The post- 
gap signal can "capture” a falling emission whose fre- 
quency lies 90 Hz below that of the post-gap signal. The 
falling emission turns around in frequency and overshoots 
toward the positive frequency direction. 

(5) The frequency difference between a gap-induced emission 
and the post-gap signal appears to be an important factor 
in determining whether or not they can mutually interact.. 


3 . Sideband Triggering 

(1) Sideband amplification and triggering are reported for 
the first time. 

(2) it has been found that the energetic electrons in the 
magnetosphere only resolve the sidebands generated by 
signals with short modulation periods. Based on this 
discovery, we have found a way to estimate the length 
of the electron interaction region. Our results show 
that it lies between 2000 and 4000 1cm. 

(3) It is found that sidebands with less than 50 Hz spacings 
mutually interact. Suppression and energy coupling among 
the sidebands are often observed. More extensive wave 
injection experiments, which will be described in detail 
in Chapter IV, reveal that 50 Hz is the order of magni- 
tude of the critical frequency range within which side- 
bands interact . 

(4) Mutual interaction between two sidebands is explained by 
the overlaps of V M ranges of electrons which the side- 
bands can organize . The electrons in the overlap can 
exchange energy with both sidebands and are responsible 
for the interaction. Therefore, 50 Hz appears to be a 
reasonable value for the corres ponding frequency range 
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of the Size of the V |( range within which electrons can 
be effectively organized by a sideband. 

(5) The sideband intensity in the interaction region is esti- 
mated to be on the order of 2.5 ~ 10 my according to a 
homogeneous model (see Appendix B) . The estimated side- 
band intensity by this simple model agrees reasonably 
well with satellite measurements Ce.g., Heyborne, 1966; 
Burtis, 1974; Inan et al, 1977] . 
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Chapter II 


VLF PULSE PROPAGATION IN THE MAGNETOSPHERE 
A . Introduction 

In this chapter, we discuss tile problem of VLP pulse distortion due 
to propagation through dispersive media . Recently, investigations on the 
interactions between whistler-mode (WM) signals and energetic electrons 
have been carried out experimentally by injecting RF pulses into the mag- 
netosphere CHelliwell and Katsufrakis , 1974; Koon et al, 1976]. Among 
the results of this work are vave growth and triggering emissions on the 
output signal as received on the ground. But the front end of a pulse 
may have been distorted by dispersion in such a way that it could appear 
to have been amplified. Without knowing how a pulse is distorted after 
traveling through the magnetosphere, we may draw wrong conclusions about 
wave-particle interaction (WPI) processes based solely on the observed 
data. Furthermore, the pulse may have been stretched enough in the front 
by the time it arrives at the interaction region to significantly alter 
the details of the WPI processes. Thus, it is essential to study the 
problem of pulse propagations in the magnetosphere. 

Pulse propagation in a dispersive medium is a well-known problem 
[Brillouin, 1960] on which numerous papers have been written. Most (if 
not all) of these studies are carried out in homogeneous media [e.g., 
Haskell and Case, 1967; Ronnang, 1968; Vidmar and Crawford, 1975; Seyler 
et al, 1972; Buckley, 1975]. As far as we know, this is the first time 
that the distortion of a VLF pulse propagating in an inhomogenous medium, 
the magnetosphere, has been calculated. 
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One of the most popular approaches in studying pulse propagation is 
to obtain an approximate analytical expression for the envelopes of the 
distorted pulses by using the saddle-point method of integration (Bril- 
louin , I960; sudden, 1966; Haskell and Case, 1967; Ronnang, 1967; Vidmar 

Ci a wf oid , 1975]. A less popular approach was adopted by Seyler et 
al C1972] to study pulse propagation numerically in a lossy but homoge- 
neous magnetoplasma. In their study, they took advantage of the avail- 
able fast Fourier transform (FIT) technique to decompose a pulse train 
into its Fourier components, then allowed each component to propagate 
through the medium separately and, finally at the points of interest, 
added all the components together by taking the inverse Fourier trans- 
form to obtain the amplitude vs time of the distorted pulse. 

Formulation of the problem of pulse propagation through a slowly 
varying medium by any of the analytical methods becomes tedious and com- 
plicated. Numerical methods using the FFT technique, on the other hand, 
require only a knowledge of the phase delays and a local modification on 
the amplitudes of those frequency components. Therefore, the numerical 
approach is more attractive than the analytical methods of solving the 
problem of pulse propagation in a slowly varying inhomogeneous medium. 

We are interested in the propagation of a whistler-mode OVM) pulse 
along a duct in the magnetosphere. By using a well-developed model for 
electron concentrations along a duct DVngerami, 1966; Park, 1971] and by 
assuming a dipole model for the geomagnetic field, the phase velocity v 
in every point along the path for each frequency component in a whistler- 
mode can be easily calculated by the dispersion relationship 


f? 


v.“ 


f(f - y 


IS 


= 1 


( 2 . 1 ) 


where f and f^ are the plasma frequency and the gyrofrequency , re- 
spectively. The phase delay can then be obtained by 



We have assumed that the spatial variation is so small that the phase 
velocity is constant within a wave length . 

In the next section, we describe the method of study in detail. In 
Section C, TO pulse propagation in a homogeneous magnetoplasma is stud- 
ied. The group velocity is reproduced. Some dispersion characteristics 
of a WM wave can be seen easily. It is observed that the distortions of 
a pulse depend on its carrier frequency. In Section D, the problem of 
pulse propagation along a duct in the magnetosphere is studied. The pulse 
shape is monitored at six locations along the path. The distortion of a 
pulse depicts the importance of the choice of the carrier frequency. We 
demonstrate that it is possible to choose a frequency at which the pulse 
envelope shows a minimum distortion at a particular location along the 
path. Discussions and conclusions are presented in Sections E and F, 
respectively. 

B. Theory 

An RF pulse has a frequency spectrum centered at its carrier fre- 
quency. The bandwidth of the spectrum, roughly speaking, equals the in- 
verse of the pulse duration. In a dispersive medium, the phase veloci- 
ties of various frequency components of the spectrum are different. As 
the pulse travels through a distance L, the phase of each frequency 
component is retarded according to its phase velocity. As a result, the 
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pulse shaDe is aH-oiwi 

- altered . It is generally (but not alWBys) true ^ ^ 
longer the propagation distance, the greater the distortion. 

We shall limit ourselves to oases in «Uch the media are linear and 
lossless. Suppose that a pulse propagates from . z = 0 to z . , and 
that the electric or magnetic field at z = 0 is 


a ^> 2 ) I z=; q - a(t,0) = q(t) sin 2f 


where 


q(t) is bound in time 
= carrier frequency 


(2.3) 


Now, let us express a (t , 


0 ) 


as a Fourier integral 


(V,°) = f A(f ,0) e 2,rift 


df 


where 


(2.4) 


A(f.O) = r a ( t|0 ) e “ 2jrift 


dt 


(2.5) 


In 0 ene3al, A(f,z) i s a complex functi 


a real and an imaginary part 


ion and can be separated into 


as 


A(f ,z) = A r (f tZ ) + i A (f ,z) 


( 2 . 6 ) 


where 


r* 00 

A r (f,z) = I a(t,z) cos 2xftdt 

CO 


(2.7) 


r m 

A i (f,z) = -J a (t ,z) sin 2*ftdt 

m/ -CQ 


( 2 . 8 ) 
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A (f ,z ) can also be expressed as an amplitude and a phase part as 


A(f,z) -A (f,z) exp iA (f,z) 
m L P h 

where 


r 2 2 -] 1 / 2 

A m (f,z) = U(f,z) + A ± (f ,z)J 


(2.9) 


( 2 . 10 ) 


V tf,z> 


tan 


-1 


A^(f ,z) 

A (f,z) 
r ' 


( 2 . 11 ) 


Unlike a pulse propagating in a homogeneous and nondispersive me- 
dium, in which the phase of its spectrum is linearly retarded (e.g., pro- 
portional to the frequencies of the components), the phase retardation 
of a pulse propagating in a homogeneous and dispersive medium defends on 
the details of the dispersion relationship. In general, the phase retar- 
dation in such a medium is not linearly proportional to the frequencies 
of the components. As a result of this nonlinear phase retardation, the 
temporal shape of the pulse becomes distorted. In addition to the non- 
linear phase retardation, there is a change in the relative amplitudes 
of the components of a pulse which propagates in an inhomogeneous and 
dispersive medium. In a slowly varying medium, this change is a local 
effect; it is discussed in Section D. 


1 * Homogeneous and Nondispersive Medium 

The terms "homogeneous" and "nondispersive" imply that the 
phase velocity is independent of space and frequency, respectively. 

Having traveled from z = 0 to z = z^, a pulse has a spectrum given 

by 
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where 


A(f,Z l ) “ A m (f ' z l ) ex P |^ A ph ^ > 2 x )] 


A m (f ' z i ) = V f ’ 0) 


( 2 . 12 ) 


(2.13) 


A ph <f,Z l > = A ph (f ' 0) ‘ ~ 


(2.14) 


Note that the phase retardation 2jtf( z /v ) is n 

v l /V p J ls l^early proportional to 

the frequency f. it can be easily shown that 


A(f,z ) - A(f,0) exp [-i2Ttf — 

\ V 


We can obtain the signal in the time domain singly by taking the 
transform of Eq. (2.15). 


(2.15) 


inverse 


a(t , 2 ^ 


•f 

•'—Of 


A(f,0) exp 


-i2itf — 
v 

P 


i2xft 


" ' * (*- %’°) 


( 2 . 16 ) 


that is, 


a(t,z,) = qft - Jl) sin 


2*f ft - -i 

• \ V P/ 


(2.17) 


Equation (2.17) indicates that the pulse has just been retarded 
by a time Z-j/v and has no distortion at all. 


2 • Dispersive Med 


mm 


The phase velocities of various frequency components are dif- 
ferent in a dispersive medium. In general, the phase velocities can he 
found through a dispersion relationship which characterizes the medium. 
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The amplitude and phase of a frequency component, after traveling from 
z = 0 to z = z^, can be written as 


A (f,z ) = A (f ,0) • m(f,z..) 
m l m l 

_ z 

= A Ph (f ' o) -I l2rtf 


(2.18) 


where m(f ,z^) is the amplitude modulation function. In a slowly vary- 
ing medium, m(f,z^) depends only on the refractive index at z^ . In a 
homogeneous medium, it is unity. 

It is noticed that the phase retardation J^ 1 2jtf (dz/v^Cf ,z)) 
is not a linear function of frequency in general . 

The spectrum of z = z^ can be expressed as 


A(f, Z;L ) = m(f, Zl ) A(f ,0) exp ^ 1 2*f v Jf )Z ) j 


(2.19) 


The corresponding signal is 


a(t, Zl ) = j m(f, Zl ) A(f ,0) exp i2itf |t - J 1 — ^ 


df (2.20) 


Equation (2.20), in general, cannot be expressed in terms of 

a(t,0) easily. By realizing that a(t f z^) is the inverse Fourier 

,zi 

transform of m(f,z^) A(f,0) exp [-i 2rtfdz/v^ (f ,z) ] , and by recol- 

lecting that A(f,0) is the spectrum of the input pulse, and that 

v (f,z) can easily be obtained from proper dispersion relations, we can 
P 

obtain a(t,z^) numerically by taking advantage of the available fast 
Fourier transform (FFT) algorithm for both homogeneous and slowly vary- 


ing media 
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We have summarized the method of the simulation of a VLP pulse 
propagating in a dispersive medium by a flow chart shown in Fig. 2.1. 

Because of the discrete nature of the operation in digital com- 
puters, analog signals must be sampled first in order to be "manipulated" 
by digital computers. In our case, both a time domain signal and the cor- 
responding spectrum in the frequency domain would have to be represented 
by a sequence of numbers properly spaced . 

Figure 2.2 shows how sampling in the time domain with a sampling 
frequency f g implies generating replicas of the spectrum in the fre- 
quency domain at the frequency intervals f [e.g., Bracewell, 1965]. 

The well-known sampling theorem stating that f must be equal to or 
greater than twice the bandwidth of the signals is sufficient to guaran- 
tee that the replicas do not overlap. We can always obtain the original 
signal in the time domain by low-pass-filtering of a properly sampled 
s igna 1 . 

Similarly, sampling the frequency spectrum of a single pulse 
corresponds, in the time domain, to generating a pulse train with a rep- 
etition rate equal to the frequency spacing between the samples in the 
spectrum. We have to work with pulse trains in our approach to the prob- 
lem. But the repetition rate of the pulses can be chosen as low as we 
want so that the pulses do not overlap. Furthermore, we shall be able 
to pick a time window, just as we picked a proper low-pass filter in the 
previous case, so as to describe a single pulse instead of a pulse train. 

Because of this discrete feature of FFT, we are physically mon- 
itoring a pulse train instead of a single pulse. The period of this pulse 
train is fixed in our simulation and equals 512 ms. The received signal 
at each monitoring station is periodic and has a period of 512 ms. 
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Therefore, by monitoring the signals at every station for only 512 ms, 
we shall have all of the information we need regarding the pulse distor- 
tion. 


C . Propagation in a Homogeneous Model 

Simulation of TO pulses propagating longitudinally is carried out 
in a model with a plasma density of 400 electrons per cc and a gyrofre- 
quency of 13 kHz for various carrier frequencies . These values are ap- 
propriate to the equatorial region near 4R in the magnetosphere. The 

E 

results have been translated down in frequency so that the carrier fre- 
quency of the input pulse is always displayed at 500 Hz. The real car- 
rier frequency may be several kHz . 

Since TO waves are circularly polarized, a three-dimensional picture 
is required to illustrate a WM wave. But we can always reconstruct a 
3-D picture by two plane-projections perpendicular to each other. Fur- 
thermore, the plan projections of a perfect circular wave differ from 
each other only in a fixed phase relationship . Therefore , only one pro- 
jection of the TO wave, as shown in Fig. 2.3, is required to illustrate 
the feature of the wave. 

According to whistler theory [Helliwell, 1965], the group velocity 
can be written as 


where 


v = 2C 
S 


i 1 / 2 ( f„ - f > 3/2 

f f 
P H 


( 2 . 21 ) 


f =s wave frequency 

f __ = gyrofrequency 
n 


fp = plasma frequency 
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The group velocity vs frequency is sketched in Fig. 2.4. The maxi- 


mum group velocity occurs at f = f R /4 and equals 


v 


\l27 


g max 


8 c f 


H 


( 2 . 22 ) 


In the model, the maximum group velocity occurs at f » 3. 25 kHz and 

has a value of 14073 km/sec An nm . 

' eC - An RF P ulse this frequency requires 

71.06 ms to travel through 1000 km. 

In the simulation , we use phase velocities to calculate the phase 
retardations for all the frequency components. We have not used the idea 
of group velocity. However, our simulation result shows that the center * 

of a RF pulse indeed propagates with the group velocity (within a mea- 
surement error of 1 percent). 

Figure 2.5 illustrates a simulation of » pulses propagating along 
the magnetic field direction. The si* panels illustrate the amplitude 
vs time records monitored at si* locations along the path. The spacing 
between the monitoring stations is 1000 km. The 50 ms pulse is injected 
into the magnetoplasma at the 10th msec at the first location (s=0km). 
It is noticed that the pulse has not been distorted very much and that 
the difference of the arrival time at the center of the pulse at two ad- 
jacent monitoring stations is about 71 ms. Because the measured delay 
and the theoretical group delay are nearly the same, we conclude that the 

wave packet is indeed traveling with the group velocity in the magneto- 
plasma . 


-- - isnrdj Jo tss s e -- - f — 


29 




Fig. 2.4. A PLOT OF GROUP VELOCITY VS FREQUENCY OF 
A TO SIGNAL PROPAGATING IN A HOMOGENEOUS MAGNETO- 
PLASMA IN WHICH f = 180 kHz AND f = 13 kHz . 

IN H 
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Figure 2.6 .ho™ the dispersion „ r 20 msec pulsea at iuimnt fM . 
quencies over a Fixed path of length (25000 ^ ^ ^ 

into the magnetoplasma at -1.526 second. The carrier frequencies are 
4.25, 3.75, 3.25, 2.75, 2.25, and 1.75 Ms. It is noticed that the pulse 
at 3.25 Ufe propagates the fastest and has the minimum distortion. We 
Shall can the frequency with minimum group delay the "nose frequency." 

In this case, the "nose frequency” is at 3.25 kHz. 

Using the leading edges appearing on the record as references, the 
arrival time could be shorter than the "group delay time” by as large as 
a few pulse widths. For example, the arrival time of the leading edge 
cf the distorted pulse in the lowest panel in Fig. 2.6 is measured to be 
-1.846 sec While the group delay time is calculated to be 1.853 sec. 
1here is a time discrepancy of about 100 msec. On the other hand, the 
arrival time of the "center" is measured to be 1.055 msec. Ihe discrep- 
ancy here is less than 1 percent. 

It is important to realize that the group velocity evaluated at the 

carrier frequency can predict the arrival time of the "center,” not the 

g edge of tile pulse. This statement applies in the frequency range 

not too close to either the gyrofrequency or zero frequency, where the 

Phase velocity equals zero, m the model we use here, the valid frequency- 

range of the carrier of a pulse longer than 10 ms Is about from 500 Hs to 
12.5 kHz . 

Other important features depicted by Fig. 2.6 are: 

(1) The distortion of a pulse depends on its carrier frequency 

the further away from the "nose frequency," the more the 

distortion. It is generally true that the distortion is 

larger be!ow the "nose frequency" than above the "nose fre- 
qucncy.” 
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1.25 kHz 


3.75 kHz 


3 . 25 kHz 


2.75 kHz 


2.25 kHz 




kHz 


Fig. 2.6. PROPAGATION OF SIX 20 ms PULSES AT 1.75, 2.25, 2.75, 3.25, 3.75, 
AND 4.25 kHz THROUGH THE HOMOGENEOUS MAGNETOPLASMA FOR 25,000 km. The 
pulses are injected at -1.526 sec. 
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C2) freq U the Carr±er 13 SbOVe the ” n ° Se fre£ J Uenc y." the low 

S frequency components arrive behind the 

h “ : ,,hen " “ 1S b81 ™ “» frequency," 

this sequence is reversed- tho uj u ^ ’ 

arrive first the 1 7 C ° mpMe « s 

til St, the low frequency components later. 

gui se P r opagation in the Magnetosphere 

In general, a wave propagatine- in nn • u 

' g an inhomogeneous dispersive medium 

ges oth m its amplitude and in its phase in 

wa . P Se ‘ In a spatially slowly 

varying medium, where the imm „ 

, vneie the WKB approximation appi ies r e _ Rl .. 

.. ^ Aes Le -g-, Budden, 1966J 

the amplitudes of e and w -p. , . ’ 

H fields are proportional to n“ 1/2 and n l/2 

respectively, so as to keep the Pnvn+- i 

„ . „ «P the Foyntxne vector |exh| constant, where 

is the refractive index rBuddcn, 1966], It is 
. . J • It is clear that the ampli- 

tudes Of E and w ^ 

depend only on the iocai refractive index not 

on the propagation path. 

In the simulation of a vtp 

S P 10 P a S a ting in the magnetosphere, 

ton xs vaixd. The refractive indices at various frequencies at a par- 
-dXar iocation of observation have been normaiiaed iocaiiy so that the 
tractive index at the carrier frequency is unity in eaicuiatin, differ- 

” tlal a " PlltUde are test in the nor- 

malization process, but the ratios between the various frequency cempo- 
nents are preserved . 

It will be sho, m later that the differential change in amplitudes 

of the frequency components in a slowlv „ 

slowly varying medium has very little 

Sffoct on the distortion of a puise. „e major contribution to the dis- 

tortion eomes from the uoniinear phase retardation in the frequency com- 
ponents . 
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The key to solve the problem Is to calculate the phase delay time 
for each frequency component. The phase velocity of any frequency com- 
ponent at a location where the gyrofrequeney and plasma frequency are 
specified is given by 


, 1/2 


v (f,z) = c 


[y a) - *| 


1/2 


f p ( Z ) 


(2.23) 


We have assumed that the refractive index is much greater than unity. 

A diffusive equilibrium (DE) model [Angerami, 1966] is used to de- 
scribe the magnetosphere in the simulation. There is empirical evidence 
in support of a DE model inside the plasmapause under normal conditions 
[Angerami , 1966; Park, 1971]. A dipole magnetic field is used to model 
the geomagnetic field. In this model, the gyrofrequeney and the plasma 
frequency at every point are given once the L value and equatorial 

density is given. L is the Mcllwain parameter [Mcllwain, 1961] 
that specifies a magnetic shell surrounding the earth. In a dipole model, 
a magnetic shell is labeled by an L when the shell’s equatorial radius’ 
equals L times of the earth radius. 

The phase velocities can be obtained for all the frequency compo- 
nents at every location along a duct through Eq. (2.23). The phase delay 
at z - z 1 can then be calculated numerically by 


V* 


'V = r 1 — 

1 Jn V 


dz 


(f ,z) 


(2.24) 


The distorted pulse at the "monitoring station" at z . , ls ob _ 
tamed by modifying the amplitude and phase of each frequency component 
according to the refractive indem and the phase delay at ^ and then 
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adding all the components together bv taking • 

S ner by taking the inverse Fourier trans- 
form via the FET. 

The inserted sketch in Fie- ? v 

r S * 2 ' 7 sho "' s «“> locations of the sin mon . 

; ° r 7 5tatiOM al0 ” g a ^ ‘ ~~ » Seated into the 

;;~ r ' “ “ 1 ~ - - «. - one pulse in every 

S. The ale length from Location 1 to To +• 

Location 6 is about 30,000 km. 

Tlie magnetosphere is modeled by the nw ^ i r 

a oy the DE-1 model (Park, 1971], and the 

equatorial density is „„ to be 400 eleotron5 ^ oc _ 

Figure 2.7 depicts the distortion of a pulse as it 
Trt „„ + . „ P 6 as it propagates from 

Location 1 to Location 6 A 10 mo •, 

P,,1Se “ 5 - 065 »• Ejected at the 

»s at location 1 . As it propagates along the duct it h 
Th(S . . . , S e duct > lts shape changes 

high frequency components appear to reach the first five „ t 
sta+ime • * five monitoring 

quenoy components lag hehind. it shonld he pointed ont that the puise 
shown in Locations I o o , . 

’ ' ’ "n 13 the same one injected at the 10th 

ms at location 1 and it arrives at Lon +■ 

Th , * 1/>0atl ° ns 5 “* S after the 512th ms. 

The pulse shorn in locations 5 and 6 is actually th 

lly the one injected into 

e magnetosphere at -502 ms from location 1. 

» la interesting to notice that the distortion is minimum at the 
equator (Location 6) at thT= 4 .- 

,. „ ““ frequency and that the travel 

time from Location 1 tn +v,„ 

the equator is shout 942 ms. According to Panes 

rSSUlt C19713 ' *• minimum group delay for a whist! 

y whistler propagating 

~ 4 m the same model magnetosphere o 

to netosphere occurs at about 5062 iiz and 

::: ™ Bat9iy i - ss ~ - — .. *** _ , inlM 
-iay occurs is called the 'We frequency,'. ^ 

" ‘ ^ th6 “ ~ - one-hop time. 

" 13 very little discrepancy between our resuit C942 ms) and Panfs 
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(940 ms). In his study, group velocities were employed to calculate the 
group delays for whistlers. In our simulation, on the other hand, the 
idea of group velocity has never been applied. It is interesting to re- 
alize that these two approaches give almost identical results for the 
delay time of a signal traveling in the magnetosphere. 

Pulses injected at other carrier frequencies, f , show different 

c 

features of distortion. Examples of 10 ms pulses at 7.075, 6.065, 4.065, 
and 3.065 Hz are illustrated in Figs. 2.8, 2.9, 2.10, and 2.11, respec- 
tively. It can be seen that the delay time from Location 1 to the equa- 
tor is indeed a minimum for the pulse at 5.065 Hz. 

When > ^ ns > the low frequency components arrive at the equator 
prior to the main wave packet and the high frequency components lag be- 
hind . There is always a location at which the distortion is minimum, as 
illustrated by the waveform at Location 5 in Fig. 2.S. The arrival se- 
quence of high and low frequency components at Location 4 is the reverse 
of that of Location 6 . 

When f c < f ng , the high frequency components arrive at the equator 
first and the low frequency components last. The further the carrier 
frequency is away from f^ , the greater the distortion. 

E . Discussion 

In addition to the "slowly varying" assumption in the simulation of 
pulses propagating along a duct in the magnetosphere, we have assumed 
that all frequency components propagate strictly parallel to a geomagnetic 
field line. Exactly longitudinal propagation seems unlikely to occur. 

In fact, various frequency components follow different snake-like paths 
inside a duct. However, in light of the theory of whistler propagation 
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[Smith, 1961], this assumption appears to be a reasonable approximation 
for wave propagation below half of the local gyrofrequeney inside a duct. 
On the other hand, exact raytracing for all frequency components inside 
a duct can, in principle, be included in the simulation if necessary. 

A duct can trap WM waves just as an optical fiber can confine opti- 
cal beams . The cross-sectional area of a duct is not uniform and is as- 
sumed to be inversely proportional to the local geomagnetic flux density. 
Therefore, wave intensity should reach a minimum on the equator since the 
largest "illumination area” occurs there. Since there is no differential 
amplitude change among the frequency components in a particular location 
due to the change in the cross section of a duct, it does not affect the 
pulse shape. Hence, it is not included in the simulation. 

Under the "slowly varying" assumption, the reflected wave energy due 
to changes in the refractive index is neglected . Thus , a wave propagat- 
ing in a slowly varying medium must have E and H fields which satisfy 
the following two physical conditions: 

(1) Constant energy flow along the propagation path. In a 
lossless medium, the Poynting vector |exh| shall be 
independent of the local refractive index, n [Budden, 

1966], 

(2) Je/H | — Z q / n, where Z q is the characteristic impedance 
of free space . 

Hence, the electric field E and the magnetic field H are propor- 
- 1/2 1/2 

tional to n and n , respectively. Note that this argument is 

not rigorous. For a more precise derivation, the reader is referred to 
Budden [1966], The condition under which the WKB approximation holds is 
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This relation is a quantitative definition of the term "slowly varying." 
It requires that the derivatives dn/dz and d 2 n/dz 2 be sufficiently 
small and that n not be too small . 

It is noticed that variations of the amplitudes of e and H are 
purely a local effect. As we mentioned before, normalization has been 
performed so that n(f Q ) equals unity in calculating the wave amplitudes 
at a particular location. The absolute wave amplitudes are lost in the 
normalization process. A true amplitude can be obtained by scaling the 

pulse amplitudes up or down by a factor of JnVT) while the shapes re- 
ma in the same . 

The contribution to pulse distortion of the differential amplitude 

changes in the WKB approximation is small in comparison with that from 

the cumulative phase retardations. A 10 ms pulse at 3.065 Hz propagating 

along a r.uct at L = 4 has been simulated for two cases: one with and 

the other without the differential amplitude changes. The result (not 

shown ) shows that the difference is negligible. This is because most 

energy (about 90 percent) of a 10 ms pulse is contained between f +100 

c 

and 100 Hz. Even at the equator, where the maximum dispersion oc- 

curs, the amplitude corrections introduced by the WKB approximation within 
this 200 Hz bandwidth are always less than 1 percent. In general, if the 
pulse is not too short (say, >5 ms) and the carrier frequency is far away 
from the local gyrofrequency and greater than, say, 500 Hz, the amplitude 
corrections can be neglected. 
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Knowing how a pulse is distorted by the dispersion in the magneto- 
sphere is essential in interpreting the observed data of many VLF wave 
experiments. An important feature of a distorted pulse is that it has 
been stretched in the front in such a way that it appears to have an 
"exponential growth" when arriving at the receiving site. For example, 
in one wave injection experiment fKoon et al, 1976], phase alternating 
signals at 6.6 kHz were injected into the magnetosphere to investigate 
the effect of WPI due to phase reversing in the triggering waves. The 
data show that the amplitudes of the 6.6 kHz signal went to zero and 
recovered with a characteristic time constant of 33 ms whenever the 
phase of the triggering wave was reversed. This time constant was found 
to be about ten times longer than the antenna current response time. 
Koon et al [1976] concluded that the recovery in amplitude was due to 
the amplifications by WPI processes in the magnetosphere. The apparent 
"growth rate" was measured to be about 260 dB/sec, which was signifi- 
cantly higher than those commonly observed [Helliwell and ICatsufrakis, 
1974]. The growth time (~33 ms) was short in comparison with that ob- 
tained from other experiments. Furthermore, there was no other evidence 
of amplification or triggering of emissions at that period. 

It was suggested by Helliwell [private communication] that the 
growth may be just another propagation phenomenon. A phase-reversing 
signal consists of several sections of RF pulses with the same carrier 
frequency. Any two successive pulses are 180° out of phase. As the 
phase-reversing signal propagates in the magnetosphere, the front and 
rear ends of each section will be distorted and stretched. Therefore, 
interference between successive pulses is expected during the transitions. 
There will be amplitude depletions during the transitions. 
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W„ have seated a phase-reversing signal at 6.6 kHx traveUng 
through « duct at h = 3.1. The equatorial plasma density is chosen to 
6, about 1000 electrons per cm 3 . These parameters are chosen to fit the 

nose frequency and minimum delay time reported by Koon et al [1976], The 

result is shown in Fig . 21 ? u 

ixg. ud. 12 . The phase is alternated at t = 0 on a 

. 100 " S PU1Se in3e ° te<1 ** - « — tion 1. As the pulse propagates 

along the duct, a "gap" is developed near the time of phase alternation. 

tVhen it reaches the senator, the gap becomes about 10 ms. We expect that 

the gap will be about 20 ms long *en it reaches the conjugate point of 

location 1. B y considering the finite antenna current response time 

(~3. 3 ms, and a finite time constant associated with the analysers, we 

speculate that the gap will appear to be 30 to 40 msec, consistent with 

Koon 1 s report 11076,. Under these circumstances, we conclude that it is 

pa B ation dxstoition, not 1VPI, which causes the amplitude depres- 

sion for about ~33 ms near the times of phase alternations. 

Another important result from these studies is that the stretching 

of a pulse is large enough bv the* 4 + 

ugn oy the time it arrxves at the interaction re- 

gion to significantly alter the details of „ processes. A pulse at 4 
KHa propagating from an altitude of 1000 hm up to the equator at h = 4 
may be stretched about 30 msec in the front, comparable to a typical 
bunching time of 20 ms (see Helliwell and Crystal [1073], Fig. 4 ). it 
is therefore necessary to either consider this distortion in the front 
of a pulse in some of the studies of WPI or build equalisers in the fu- 
ture to compensate for sue, distortion. »or example, the experiment of 
Phase alternations can conceivably be carried out in the future by "pre- 
processing” the injected signal so that, when a phase-reversing signal 
reaches the interaction region, the phase of the signal will appear as 
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being reversed instantaneously. Figure 2.13 illustrates a pre-processed 
signal propagating along the duct at L ~ 3.1. As the signal readies the 
equator, the signal becomes 'clean.” By the same principle, post— proces- 
sing of the received signal can be implemented to compensate for the dis- 
tortion due to propagation from the interaction region to the receiving 
site . 

In the experiments designed to study wave-wave interaction (WWI) in 
the magnetosphere, frequency shift -keyed (FSK) signals are injected into 
the magnetosphere [Chang and Helliwell, 1B77J. A FSKed signal consists 
of RF pulses alternating between two carrier frequencies at regular in- 
tervals whidi are usually in the multiples of 10 ms. The sidebands gen- 
erated by the FSKed signals are constant frequency signals and last as 
long as the FSK signals continue, usually one second or longer. In other 
words, each sideband is an RF pulse which lasts one second or more in- 
stead of on the order of 10 ms. The front and rear ends of sidebands 
are not of our concern. We are using the middle portions of the pulses 
to investigate WWI . Thus , the distortion of pulses as a result of prop- 
agation will not jeopardize our conclusions in the studies of sideband 
mutual interactions. 

F. Conclusion 

Among the results of the study of pulse propagation in the magneto- 
sphere are the following: 

(1) By properly choosing the carrier frequency of injected 
pulses, the propagation effects can be minimized. The 
best frequency to choose is the "nose frequency" (as 
measured at the equator) . 
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(2) in designing certain wave-in jeetion experiments in 

the study of WPI processes in the future, such as the 
phase-alternation experiments, we should include pre- 
and post-processing filters in the transmitting and 
the receiving facilities, respectively, to remove the 
distortions due to propagation in such a highly dis- 
persive medium. 
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Chapter III 


GAP INDUCED EMISSIONS 


A. Introduction 

Recently, many VLF wave injection experiments have been performed 
to investigate the details of VLF wave-particle interaction (WPI) pro- 
cesses in the magnetosphere. As a result of this interaction, coherent 
VLF waves may grow exponentially with time and trigger a variety of emis- 
sions [Helliwell and Katsufrakis, 19743. it is generally believed that 
the amplification and triggering mechanisms are due to cyclotron reso- 
nance bt ,veen energetic electrons and VLF waves. Beyond that, there is 
no general agreement. 

It has been suggested that coherent wave amplification is the result 
of the phase bunching of electrons [e.g., Helliwell, 1967]. The TO waves 
phase-bunch the electrons, producing transverse currents that cause the 
waves to grow. In Appendix B, we shall discuss the mechanisms of VLF 
wave-particle interaction in the magnetosphere to give the background 
and rationale for many VLF wave injection experiments related to this 
report . 

In Section B, we present an unejqDected and interesting result from 
one such experiment designed to study how the wave growth is affected by 
altering the phase of the triggering signal. The experiments have been 
conducted between Siple, Antarctica, and its conjugate point near Rober- 
val, Quebec [Helliwell and Katsufrakis, 19743. Because of inherent lim- 
itations of the Siple transmitter, a phase alternation in a VLF signal 
is achieved by inserting a 10 msec pulse offset in frequency from the 
main signal by several hundred hertz . The frequency offset is so large 
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that the 10 msec offset pulse does not interact significantly with the 
electrons in resonance with the main signal. The main signal can there- 
fore be viewed as having been interrupted by a 10 msec gap. The portion 
of the signal behind the gap is either in phase or in antiphase with the 
one preceding the gap. 

The data show that the post-gap growth rates and the saturation 
levels of the main signals are the same regardless of whether the phase 
of the triggering waves have been shifted by zero or rt-radians during 
the 10 msec gaps. Furthermore, a rising emission often is induced by 
the 10 msec gap. This gap-induced emission closely resembles the tail- 
end emissions that frequently appear at the ends of constant frequency 
pulses. These gap-induced emissions may develop into fully independent 
rising emissions or they may be either suppressed or entrained by the 
main signals. 

The key to an explanation of the phenomenon of gap— induced emis- 
sions is based on the hypothesis that the radiation from the wave orga- 
nized electrons can switch from a forced mode to a natural mode at the 
end of the triggering wave. Triggered emissions are said to be in a 
natural mode because they are not driven by external sources . The ampli- 
fied signals on (or near) the frequencies of triggering waves are in a 
forced mode since they are driven by the triggering waves. The detailed 
mechanism of the development of a natural mode emission is dealt with 
separately [Helliwell, 197S] . a qualitative explanation is presented in 
Section C. 

A pi eliminary result from a new wave injection experiment (gap— trig- 
gering experiment) will be presented in Section D. This experiment has 
been designed to investigate how large a gap shall be to allow a falling 
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emission to become fully developed wlthln the gap . 0n8 Qf ^ 

confirms an earlier observation that falling emissions always start with 
a small rise in frequency femes, 1S74J. A new feature of WI „ iden _ 

tified. A post-gap signal can capture a gap-induced falling emission 

whose frequency is 90 Hz below the post—au „ • 

S P gnal, causing the emission 

to turn around toward the positive frequency direction. 

in Section E, the results are solarised and the recommendations 
for future studies on this subject are presented. 


B ‘ Examples of Gap-Induced Emissions 

It is assumed that TO wave amplif ication and triggering in the 

magnetosphere are due to nonlinear phase-bunching of energetic electrons 

by VLF waves. The "nhasoc" ~-p n 

of electrons refer to the angles between the 

transverse velocity components of the electrons and the wave magnetic 
Held, it was thought that the phase-hunched currents might be diminished 

^ Pertm ' bi ” g thS PhaSM ° f ^ — s. A transmitter program 

“ ~ d *° S * Udy and the saturation ievel of a 

wave would be affected simply by reversing the phase of the triggering 

e idea is shown in Fig. 3.1. Electrons are randomly distributed 
m phase Initially, The wave organizes ths eleotrons by y XB 

The electrons are phase hunched after a period of time. At" tit, we 

reverse the phase of the waves. Whether the 1 

> nether the electrons would be debunched 

by this reversed-phase wave is the question. 

Because of inherent limitations of the Siple transmitter, a 10 ms 
pulse offset from the main signal at f p by a few hundred hurts must be 
inserted in order to aobieve a phase change in the main signal. The 
amount of phase change in the main signal is determined by 
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a . Electrons initially 

are uniformly dis- b. After a period of time, 

tributed in phase. electrons are phase- 

There is no net organized to form a 

transverse current . transverse current J . 

x 



c. The phase of the sig- 
nal is reversed. The 
phase-organized elec- 
trons may suffer 
opposite forces and 
undergo a debunch ing 
process . 


/ 

/ 

/ 



B w 


d . After the electrons 
become randomized in 
phase, the post-tran- 
sition signal starts 
afresh to organize 
these electrons . 


Fig. 3.1. DEMONSTRATION OF THE IDEA OF HOW A PHASE REVERSING 
SIGNAL AFFECTS THE PHASE-BUNCHED CURRENT. Short arrows rep- 
resent Vj^'s of the resonant electrons. It is assumed that 
v ± is constant and that electrons initially are in exact 
resonance with the wave . 
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(3.1) 


= 2jtAfAt 

A schematic illustration of how to achieve a phase change in the 
roain signal is shown in Fig. in this experiment, * ia fixed at 

10 " S - SOtti " e " *° * 300 “ “ 5 ° »■ .1". a phase shift of * „ and 
-7n, respectively. Bins, the in-phase and the anti-phase conditions 

are achieved, for simplicity, we shall call the waves in the in-phase 
condition O-waves and those in the anti-phase condition t-waves. 

At the time the experiment was designed, it was believed that the 

10 “ " P bet °° “ — '-e growth processes and that 

the main effect would come from the phase reversing. T he data show Just 

the reverse. Bie growth rates and the saturation levels of the signals 

after the interruptions are the same regardless of whether the triggering 

wave is a a-wave or an C-»ave. in addition, rising emissions often are 

y the gaps. These gap-induced emissions may develop into fully 

independent rising emissions or they may he either « 

i y ne either suppressed or en- 
trained by the main signals. 

TWO exampies of the gap-induced emissions are illustrated In Plg 
3-3. Rising emissions are developed from the 10 ms gaps. B,e fluency 
and time resolutions of these records are about 20 Ha and 60 msec, re- 
spectively. The dynamic speotrums in the two lower panels, showing the 
seme data as the ones on the upper panel, have been expanded in fluency 
by a factor of 2 and 5, respectively. Bie two triggering .eves at 5 RHs 

* WaVSS ‘ B ° th laSt f ° r ° ne SeCond and are interrupted at the 400th 

' SM by “ 10 “ gaP ' « th. signals are indicated by the 

one-second long horisontal bars Just above the time scale. Bie time of 

interruptions are indicated by the email vertioai bars. Bie delay time 

of the gap with respect to the beginning e, the signal is called T 

a ' 
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Pig. 3.2. DEMONSTRATION OP HOW A PHASE REVERSAL IS ACHIEVED. 
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frequency scale by factors of two and five respectively' TothT^T ^ 
triggering signals at 5 kHz are interrunted hv l in ^ * B h f the ° ne_s 

£ r ratl ° n r SlE,,alS “* indicated beneath* the "Z IZZTZZollZ 

?z::i °ztz: lndica r the time ° f 

and sr 

dt ~ 1 kHz/sec. ' ’ . m turn , induce rising emissions with df/ 
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There were several ducts in the magnetosphere in this period as 
evidenced by multiple traces of whistlers on the Roberval data (not 
shown) . It is believed that each whistler trace observed on the ground 
corresponds to a physical duct in the magnetosphere [e.g,, Carpenter, 
1966; Park, 1970], Only two whistler traces were observable above 4.S 
kHz. The two-hop delays were measured to be ~2.10 and ~2.62 sec at 5 
kHz for these two whistlers. Moreover, the second whistler appears to 
be much intenser in amplitude than the first one. (The L values of these 
ducts were found by a standard curve-fitting method Ce.g., Smith and 
Carpenter, 1961; Ho and Bernard, 1973] to be at L = 3 .3 and L = 3.5.) 

The one-s pulses at 5 kHz were transmitted on the seconds from 
Siple. The one-hop delay of the signals (referring to the gaps and the 
ends of the pulses) was measured to be ~ 1.3 sec, in good agreement with 
the delay time of the second whistler (at L ~ 3.5). Therefore, we con- 
clude that the received signals can be viewed as single-path signals. 

The first example shows that a rising emission is developed near 
the time of the 10 msec gap. The slope of the gap-induced emission is 
about 8 kHz/sec. As the emission develops, it appears to feed energy to 
the 85th and 86th power line harmonics (PLH’s) at 5.10 and 5.16 kHz , 
respectively. The PLH's then trigger emissions with df/dt « 1 kHz/sec. 
The gap-induced rising emission with df/dt «8 kHz/sec is quenched near 
5.25 kHz . 

The PLH's are believed to be signals radiated by the North American 
power distribution systems leaking into the magnetosphere CHelliwell et 
al, 1975]. 

The second example shows almost the same features of gap-induced 
emissions as the first example. The gap-induced emission with df/dt ~ 8 
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kHz/sec is quenched near 5.3 kHz. Only one PLH, at 5.1 kHz, has been 
activated by the gap-induced emission and triggers a rising emission with 
slope of ~1 kHz/sec. Near the time of triggering, the main signal sud- 
denly decreases in magnitude. As the rising emission drifts away from 
the main signal, the main signal starts to grow again. 

The amplitude behavior of the main signal before and after the gap 

is illustrated by four more examples, shown in Fig. 3.4. Here, t is 

d 

200 or 400 ms, as indicated in the middle panel by a vertical bar below 
the ”0” or "it" symbol. 

It has been determined through whistler studies that there were 
several paths available at that time. The dominant one was at L ~ 4.6. 

The signal amplitude measured in a narrow-band filter (340 Hz) cen- 
tered at 4.5 kHz shows an initial growth rate of (75 ± 10) dB/sec in all 
four examples and a growth time of roughly 250 ms. It is seen that, 

when r =200 ms, the gap occurs before saturation while, for t =400 
a d 

ms, the gap occurs after saturation. There are very pronounced amplitude 
drops associated with the gaps in the main signals. The growth rates 
after the gap are the same as the initial growth rate (within a limit of 
observation error), regardless of whether the signal is an 0-wave or a 
rt-wave . The 0 and it conditions do not seem to affect either the 
growth rate or the saturation levels after the gap. 

The spectrogram in the upper panel of Fig. 3.4 shows that the emis- 
sions triggered by the gaps are not fully developed. They may have been 
suppressed by the main signal or by some other process. The evidence 
suggests that emissions are, in fact, suppressed by the main signal. 

This suppression effect is demonstrated by the records on Fig. 3,5. 
The two examples of gap-induced emissions in the middle panel occurred a 
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Fig. 3.4. EXAMPLES OF GAP-INDUCED EMISSIONS. Upper panel shows the dynamic spectrum. Lower 
panel shows amplitude (log scale) measured in a 340 Hz band centered at 4.5 kHz. The in- 
terruptions occur at the 200th or 400th msec as indicated by the vertical bars in the mid- 
dle panel. Shortly after each gap, the amplitude drops approximately to the noise level 
and then grows again. The initial and post-gap growth rates are almost identical regard- 
less of whether the triggering wave is in an 0- or a rt-condition . The saturation level 
is also independent of the 0 or « condition. The growth time is about 250 msec. 
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few seconds earlier than the ones displayed in Fig. 3.4. Here, 
equals 200 and 400 ms, respectively, as indicated by the vertical bars 
just above the panel. Hie upper and lower panels show the calibration 
waves prior to and following the examples shown in the middle panel 
(within 30 sec). Hie calibration waves, as indicated by horizontal bars 
above the panels , are constant frequency pulses of length 200 and 400 
ms, exactly the same as the pre-gap sections of the waves in the middle 
pane 1 . 

It is clear that both the 200 and 400 ms calibration waves trigger 
"band-limited-impulse" (BLI) type emissions [Helliwell, 197S] , followed 
by fully developed falling and/or rising emissions. But the pre-gap 
sections (200 or 400 ms long) of the gap-triggering waves trigger BLI- 
type emissions only. There are no fully developed falling or rising 
emissions following the BLI's, suggesting that the emissions are either 
suppressed or entrained by the post-gap section of the main signals. 
Remembering that the post-gap growth rate of the main signals is the same 
as the initial growth rate in all examples shown in Fig. 3.4, we rule out 
the possibility of entrainment. Therefore, we conclude that the emis- 
sions have been suppressed by the main signals and that the main signal 
following the gap has not been significantly altered by the emissions 
in this particular case , 

Table 3.1 lists the averaged initial growth rate, post-gap growth 
rate , and the saturation level over 52 examples taken on April 28 and 
May 13, 1975. The result shows that the saturation levels are practi- 
cally the same and that the initial growth rate is identical to the 
post-gap growth rate, regardless of whether the wave is in a 0- or a 
condition . 
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Table 3.1 

AVERAGED INITIAL GROWTH RATE, POST-GAP GROWTH RATE, 
AND SATURATION LEVEL OVER 52 EXAMPLES TAKEN ON 
APRIL 2S AND MAY 13, 1975, FOR BOTH O-WAVES AND jt-WAVES 


T d 

200 ms 

400 ms 

Avei’age 

Wave type 

0 

it 

0 

it 

Initial growth 
rate (dB/s) 

79 

(26)* 

79 

(23) 

76 

(27) 

7S 

(23) 

7S 

Growth rate 
after gap (dB/s) 

74 

(27) 

78 

(26) 

74 

(30) 

77 

(23) 

76 

Saturation** 
level (dB) 

15 

(2.8) 

14 

(3.3) 

16 

(2.S) 

16 

(3.2) 

15 


(M) indicates the standard deviation. 
Above noise level . 
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The data on gap-induced emissions show many instances of rising 
emissions. So far, no fully developed falling emission has been observed 
within the gap in more than a hundred cases we have studied. This indi- 
cates that falling emissions require a longer time to develop. This is 
understandable since all tail-end falling emissions first show a slight 
rise before falling [Stiles and Helliwell, 1975], By the end of the 
gap, the frequency of the falling emission is still within the range for 
which interaction between the post-gap signal and the emission can occur. 
Thus, the falling emissions are suppressed by the post-gap signals. 

C. An Interpretation of Gap-Induced Emissions 

The key to an explanation of gap-induced emissions is based on the 
hypotheses that the radiation from the wave -organized electrons can 
switch from a forced mode to a natural mode at the end of a triggering 
wave. As the wave propagates through the equatorial interaction region 
in the magnetosphere , it phase-bunches electrons , producing transverse 
currents that radiate a new field, causing the wave to grow. The fre- 
quency of the new radiation field usually remains within a few hertz of 
the triggering wave [Stiles and Helliwell, 1975]. The new outgoing 
radiation organizes incoming electrons, which then radiates new fields 
again and the process continues. As noted by Helliwell [197S] , this 
"wave-particle system" resembles a resonant circuit driven by signals 
which may not be at the natural frequency of the circuit. The resonant 
circuit is in a driven mode or forced mode. When the external signal 
disappears, the resonant circuit switches to its natural resonance fre- 
quency and is then operating in a natural mode. 
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As a gap-triggering wave arrives at the equator, the gap has ohanged 
because ol the fact that the wave has been traveling through a highly 
dispersive medium. Computer simulation results (not shown) indicate 

that the gap is expanded to 20 to 25 „ 5 for the cases discussed in the 
previous section. 


When the electrons organized by the pre-gap section of a wave cross 
the gap (this transit requires about 10 ms because their parallel veloc- 
ity makes the gap appear shorter than the lab frame value of about 20 
ms), there are several possible responses for the "wave-particle system." 
First, when the gap is very small compared to the debunching time, the 
coherence of the organized electrons tends to be preserved during the 
gap. Assuming no inhomogeneity in the interaction region, the post-gap 
section of the wave is able to control these previously bunched electrons 
if the triggering wave is an 0-wave. The signal then continues to grow 
when the interruption occurs in the initial growing phase. On the other 
hand, if the wave is a rt-wave, the post-gap section has to reverse the 
Phases of currents, if that can be done, before the wave grows again. 

We expect to see a systematic difference in growth behavior after the 
gaps for these two cases . 

An electron motion in an 0-wave and that in a rc-wave is shown in 
Figs. 3.6 and 3.7, respectively. On the left of both figures, an elec- 
tron trajectory in the phase plane (V ( ,i) prior to, during, and fol- 
lowing the gap is shown in panels (A), (B), and (C), respectively. The 
corresponding time domain motion is shown on the right. Prior to the 
gap, the electron is trapped by the wave and executes a pendulum type 
of motion as indicated by 1, 2, and 3 in both the phase plane and the 

time domain. During the gap , V(| of the electron is constant because 
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Fig. 3.7. ELECTRON MOTION IN A it WAVE. (A) Before the gap, (B) during 
the gap, and (C) after the gap. The electron becomes "untrapped" after 
the gap . 


there is „„ wave that ea„ effectively perturb the electron. The phase 
angle * changes at a constant rate KCV, -V ). In the phase plane, the 
electron traces a line parallel to *-a*is while, in the tine domain, it 
rotates with a constant angular rate of KlV^-y/a* cycles per second. 
At the end of the gap, the electron reaches Location 4. After the gap, 
the electron can be trapped in the O-wave but not in the t-wave . In the 
O-wave case, the electron has been shifted to a different contour closer 
to the origin in the phase plane. The corresponding tine domain notion 
has been changed too. The electron still swings back and forth about the 
wave field but with less excursion in * .„d with a smaller oscillation 

period, in the t-wave case, the electron has been moved to an unclosed 
contour in the phase plane, as shown in Pig. 3.7. it beco.es untrapped, 
in the time domain, the electron rotates about the wave field. Thus, this 
particular electron remains trapped after the gap for the O-wnve case and 
becomes untrapped for the jt-wave case. 

Not all the trapped electrons become untrapped after the gap it. a 
*-«uve. omy a portion of them do. Figure 3.S illustrates the temporal 
development of the trapping region during the gap. me inhomogeneity 
forces have been neglected for the time being. Panels (a), (b), and (e) 
show the cases In which the gap is 0, 10, and 20 msec, respectively. 
O-wave and tr-wave conditions are indicated by the V and 'V symbols. 

The regions enclosed by the broken lines are the deformed trapping regions 
developed during the gap. Assuming the wave intensity is the same for 
both the pre- and the post-gap sections, the trapping regions for the 
post-gap signal are indicated by the solid lines. The overlap of the 
deformed and the new trapping regions is shorn by the shaded areas in 
which electrons trapped by the pre-gap signals are trapped again by the 
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post-gap signals. Those electrons outside the shaded area become un- 
trapped and are assumed to make no contribution to the growth of the post- 
gap signals. It must be noted that using trapped electrons to calculate 
the transverse current may only give us a rough estimation. However, by 
examining how many percent of the trapped electrons have been lost during 
the gaps, it is possible to determine whether the post-gap wave will take 
advantage of the previously organized electrons or start afresh to orga- 
nize a new set of electrons. In this illustration, the wave intensity 
and particle pitch angle are set up so that the corresponding frequency 
range of the trapping width in V ( equals 50 Hz. it is observed that 
the 0-wave has more control over the previously organized electrons than 
the jt-wave for a small gap (less than 10 ms). 

The inhomogeneity force can be taken into account by a model of a 
constant external torque, as discussed in Appendix B. This torque deforms 
the trapping range, reducing the size of the range and moving the stable 
point f = 0 to other values of f In addition, the V (| changes dur- 
ing the gap because of the adiabatic force. We have used a simple model 
for the inhomogeneity by considering a constant V ( change during the 
gap for all the trapped electrons, neglecting the deformation effects in- 
troduced by the external torque. But the deformation by the differential 
t-drift due to V ( spreading is included. The result is shown in Fig . 

3.9. The trapping width in V (| is chosen so that the corresponding range 
in frequency is 50 Hz . The inhomogeneity has changed the electron V 's. 
Usxng typical parameters at L = 4 and assuming electron pitch angle at 
30 , we can obtain a V |( change that corresponds to a 15 Hz change in 
frequency at a distance 500 km away from the equator within 5 msec. It 
is clear that about 50 percent of the trapped "area" becomes "untrapped" 
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Pig. 3.9. DEFORMATION OP THE TRAPPING DOMAIN DURING 
A 10 msec GAP IN AN INHOMOGENEOUS MODEL. The tran- 
sit time of an electron to cross the gap requires 
about 5 ms because the electron and the wave travel 
in opposite directions. 
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frequencies of the organized electrons are shifted from the main signal 
frequency at a rate corresponding to the local inhomogeneity. After the 
gap, there are two signals: an emission and the main signal. When the 

frequency difference between the two signals is large, they behave inde- 
pendently and the emission can become fully developed. At the same time, 
the main signal organizes its own electrons from the beginning, regard- 
less of whether the signal is an O-wave or a it -wave. On the other hand, 
when the frequency difference is small, the emission may interact with 
the main signal. In this situation, the emission may be entrained or 
suppressed by the main wave. The emission may also suppress growth of 
the main signal. 

The data discussed in Section B show that a 10 ms in a triggering 
wave can induce emissions in favor of the third possibility. The remain- 
ing questions are: 

(1) What are the criteria for a self-excited oscillation in 
tlie wave -particle system? 

(2) How long does it take for radiation from this oscillat- 
ing system to switch from a forced mode to a natural 
mode at the end of an external wave? 

(3) For two waves to interact, are there conditions between 
the relative phase, wave frequencies, and wave ampli- 
tudes to produce different effects such as entrainments 
and suppressions? 

These questions have not yet been answered satisfactorily. However, 
the gap triggering data and related experiments do give us some insights 
on the questions. 
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It should be recalled that short individual pulses « 50 ms) do not 
usually trigger emissions, while the long pulses (>100 ms) do Dlelliwell, 
1965; Stiles and Helliwell , 1977 ; Helliwell and Katsuf rakis , 1974 ; McPher- 
son et al, 1974]. This suggests that the wave amplitudes must grow 
beyond a threshold level in order to generate self-sustaining emissions . 
The question of how the threshold levels depend on input wave intensities 
is still under investigation. 

In a simple lumped resonant circuit, the oscillation switches from 
the forced mode to a natural mode "instantaneously" when the external ex- 
citation is removed. On the other hand, the wave particle system in the 
magnetos phere may not be able to switch from a forced mode to a natural 
mode instantly. One of the reasons is the dispersion distortion on the 
ends of signals. A VLF pulse may appear to have a "tail" when arriving 
at the interaction region near the equator. The "end" is not a sudden 
termination. Since a 10 ms gap in the triggering wave can induce emis- 
sions, we conclude that the radiations from the organized electrons can 
switch from a forced mode to a natural mode in less than 10 ms. 

After the gap, the main signal and the triggered emission interact 
with incoming electrons at the same time. If the frequency difference 
is large, these two signals will organize their own groups of electrons 
independently. They do not interact with one another. On the other hand, 
if the frequency difference is small, interaction between the gap-induced 
emission and the post-gap signal may occur. 

The detailed mechanisms of this kind of TO I are not yet clearly 
understood. It is, however, believed that both the inhomogeneity of the 
geomagnetic field and the overlap of the coherence bandwidth (see the 
next chapter) are responsible. In the next chapter, we shall show that 
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two forced-mode signals with frequency difference less than 50 Hz tend 
to interact with each other. Coupling, suppression, and entrainment are 
often observed in such signals . The interaction is attributed mainly to 
the overlap of perturbed V ranges in the electron distribution func- 
tion [Chang and He lliwell, 1977], Electrons in the range of overlap 
carry energy between these signals. The same mechanism may happen m 
FNWWI. This hypothesis needs further investigation. On the other hand, 
Helliwell [private communication, 1978] suggests that the inhomogeneity 
is the main source of wave couplings . 

As to the question of whether the frequency slopes and the relative 
phase between two interacting waves are important to the processes of 
entrainment, suppression, or coupling, the available gap-triggering data 
are not able to provide a proper answer. More refined, controlled exper- 
iments or computer simulations are required to solve this problem. One 
of the refined experiments, Gap Triggering Experiment, will be discussed 
in the next section. 

Xn summary, the phenomenon of gap-induced emissions is explained by 
the hypothesis that radiation from the "wave-particle system" switch from 
the forced-mode to a natural mode at the end of a triggering wave. Natu- 
ral mode emissions are developed in a 10 ms gap. The wave intensity is 
strong enough to maintain a self-sustaining emission. The frequency of 
the emission during the gap changes according to the local inhomogeneity 
in the interaction region. The frequency difference between the emission 
and the post-gap section of the main signal may be large after the gap. 
Tne emission may then have developed into an independent rising or fal- 
ling tone, and the post-gap section of the main signal starts afresh in 
organizing new sets of electrons. The 0- and ic-condition is immaterial 
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as fr.v as the growth of the post-gap section of the signal is concerned. 
On the other hand, when the frequency difference is small, interaction 
between the gap-induced emission and the post-gap signal occurs. The 
frequency difference between the emission and the main signal after the 
gap is the main criterion in determining whether mutual interactions 
occur or not . 

D . GATH Transmitter Program 

The gap-triggering data show many instances of inducing rising emis- 
sions from the 10 ms gaps. So far, no fully-developed falling emission 
has been observed within the gap, indicating falling emissions require a 
larger gap to develop. It is believed that a fully developed falling 
emission might have been observed, had the triggering signal had a larger 
gap. It is therefore beneficial to increase the length of the gap to de- 
termine the critical length beyond which a fully-developed falling emis- 
sion can occur. We expect that waves with gaps less than the critical 
size will produce mi's. This program also provides an opportunity to 
actually scan the development of falling emissions. These are the 
reasons for designing the so-called GAP-TRIGGERING (GATR) PROGRAM. 

In all the examples shown in Fig. 3.10, the 750 msec triggering 

waves are interrupted by a gap at the 250th msec. The gap size, T , 

g 

has five different values: 10, 40, 70, 100, and 130 msec, and all the 

waves are in the 0-condition (no phase reversal). The carrier frequency 
is 3.S5 kHz. From the multiple traces of the ramps appearing on the 
records (not shown) , it is found that there were two dominant paths at 
this period. They were very close to each other. The difference of the 
time delay for these two paths at 3.85 kHz is about 50 ± 10 ms (refering 
to the centers of the traces ) . 
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Fig. 3.10. THE FEATURES OF GAP- INDUCED EMISSIONS VS 
VARIOUS GAP SIZES. 


It is observed that , as the gap increases from 10 to 70 ms , 

all the gap-induced emissions are either suppressed or entrained by the 
post-gap signals. When T g is 100 ms or more, a fully-developed fal- 
ling emission is observed from one of the two paths. On the other path, 
path 2, even when the gap becomes as large as 130 ms, the emission can 
still be 'captured" by the pcst-gap section of the main wave. This is 
because the slope df/dt of the emission is small on path 2. After the 
130 ms gap, the frequency difference between the emission and the post- 
gap signal is small enough to permit the WWI to occur. This "capture" 
feature can also be seen for the signal on path 1. When the gap equals 
70 ms, the gap-induced falling emission has been developed. It rises 
in frequency initially and then turns around in frequency as a typical 
tail-end falling emission. At the end of the gap, the emission frequency 
is 90 ± 10 Hz less than that of the main signal. The post-gap signal 
captures the falling emission, causing the emission to overshoot toward 
the positive frequency direction. This "capture" feature of WWI is id- 
entified for the first time through this program. 

Explanations of this "capture" feature of WWI are still being in- 
vestigated. It may be an important clue to understanding the mechanisms 
of WWI in the magnetosphere . 

E • Summary and Recommendations 

We have found that a 10 ms gap in a triggering wave can induce 
emissions. The relative phase between the pre- and post-gap sections is 
immaterial as far as the growth processes of the post-gap sections are 
concerned. A qualitative explanation has been found by postulating that 
the radiation from the organized electrons can switch from a forced-mode 
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to a natural mode at the end of a wave train. The 0- and ^-conditions 
are not important to the growth of post-gap sections of signals. On most 
occasions, the gap-induced emissions interact with the post-gap signals. 
This wave-wave interaction involves one natural mode emission and a 
forced-mode signal and is termed force mode/natural mode wave-wave in- 
teraction (FNWWI). For cases shown in Figs. 3.4 and 3.5, we have shown 
that the growth of the post-gap signal seems not to be affected , at least 
not noticeably, and that the development of the natural emission is 
quenched by this interaction, regardless of whether the wave is in an 
0- or a jT-condition . 

A z’efined transmitter program based on the observations of the 
gap-induced emissions has been designed to actually "scan" the develop- 
ment of triggered emissions . One of the results confirms an earlier 
observation that triggered emissions always start with a small rise in 
frequency. It has also been found that falling emissions with a negative 
frequency offset as large as 90 Hz from the triggering wave can still be 
"captured" by the post-gap signal. This is the first time the "capture" 
feature of FNWWI have been identified. This refined transmitter program 
can produce FNWWI' s in a somewhat predictable way. For small gaps, sup- 
pressions are usually observed . For large gaps (>100 ms), falling emis- 
sions can be fully developed. For medium size gaps (~70 ms), one form 
of entrainment, the "capture" feature, is often observed. It seems feas- 
ible to use this program to generate various features of FNWWI 's . A re- 
fined version could be employed to investigate the importance of the rel- 
ative phase, the frequency difference, and the slope df/dt of the 
emission in a FNWWI by controlling the phase and the frequency of the 
post-gap section of the triggering wave. 
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Chapter IV 


SIDEBAND TRIGGERING 

A. Introduction 

In this chapter, we shall discuss another kind of wave-wave inter- 
action (WWI): sideband mutual interaction (SMI), it was discovered 

unexpectedly during a wave injection experiment for studies on the re- 
lation of signal growth to phase reversals in triggering signals. The 
triggering signal can also be categorized as an O-wave or a ir-wave, as 
before. The rationale for this kind of wave injection experiments has 
been discussed in the preceding chapter. In this particular experiment, 
which will be described in detail in Section B, the phase of a one-sec- 
ond signal is altered regularly at periods that are multiples of 10msec. 
At the time of designing the experiment, we did not consider the discrete 
sidebands at various frequencies that were generated by the regular phase 
modulation of the signals. We were expecting to see systematic differ- 
ences between the growth of 0-waves and that of *-waves . The data show 
that 0- and rt-cond it ions are immaterial to wave growth and triggering. 
Furthermore, the period of phase modulation in the signal appears to be 
closely related to the occurrence of triggered emissions. The longer the 
period, the less the occurrence of triggered emissions when the period 
is less than 100-200 msec. Signals with modulation periods longer than 
200 ms trigger emissions similar to those triggered by one-second pulses. 
THese observations lead to the discovery of sideband triggering. 

In Section c, we shall discuss how the length of the interaction 
region can be estimated from our observations. This estimation is n 
measure of electron "memory time" which is assumed to be the same as 
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the transit time for electrons to pass through the interaction re- 
gion. 

In Section D, we shall show that the emissions triggered by indi- 
vidual sideband components are similar to those triggered by constant 
frequency signals. This phase modulation technique results in the in- 
jection of multiple waves (sidebands) into the magnetosphere. The fre- 
quency spacing between these sidebands have been changed by varying the 
modulation period . Therefore, this experiment becomes the first con- 
trolled experiment on SMI in the magnetosphere. We shall also show two 
features of energy coupling between two sidebands at different frequencies. 

In Section E, we shall discuss the mechanisms of TO I in both time 
and frequency domains. Then, a simple model of coherence bandwidth is 
used to estimate the wave intensity in the interaction region. 

In Section F, we shall show results from a new transmitter program 
which was developed recently based on the studies of the first experi- 
ment on SMI. We shall show that, for most cases, the coherence bandwidth 
is about 50 Hz. The corresponding wave intensity in the interaction re- 
gion is about 2 . 5 ~ 10 my, in reasonable agreement with satellite measure- 
ments tHeyborne , 1966 ; Burt is , 1974; Inan et al, 1977]. More examples 
related to the coherence bandwidth from various transmitter programs are 
shown in Appendix D. 

The conclusions are presented in Section G. 

B . A Transmitter Program for Generation of Sidebands 

The special transmitter program described below exhibits several 
features, including phase shifts, shifts in frequency, and variations 
in pulse length. The phase shift of a signal at a given frequency is 
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The 


accomplished by shifting the frequency by Af for an interval x 
phase shift £0^ equals 2xAff s . In this experiment, Af ranges be- 
tween 100 to 150 Hz and t s ranges from the minimum achievable length 

10 to 100 msec. The phase is shifted at a regular interval of 2 T at 

s 

a given frequency. Figure 4.1 schematically sho.ws part of the waveform 

and a frequency-time (f-t) display of a one-s wave with r =10 msec. 

s 

The total program lasts for 30 seconds. It consists of 30 one-s 

waves which are grouped into five sets according to the value of r . T 

s s 

has the values of 10, 20, 30, 50, and 100 msec. Figure 4.2 illustrates 
two of the five sets with = 50 and 100 msec. Each set has six waves 
which alternate between two frequency bands separated by 500 Hz. The 500 
Hz separation applies to the lower frequency of the upper band and the 
upper frequency of the lower band . Hence , the offset frequencies Af are 
positive in the upper band and negative in the lower band. In each set, 
jAf| exhibits three different values: 0, 100, and 150 Hz, except when 

T s = 20 and 100 msec. In these cases, (Af j are 0, 100, and 125 Hz. 
lVhen Af = 0 Hz, the signal is a one-s constant frequency wave and is 
called the calibration wave or C-wave. When Af = 100 Hz, the 0-wave 
or in-phase condition is produced and, when Af = 125 or 150 Hz, the 
x-wave or antiphase condition is achieved. The transmission sequence 
for the upper frequency band in a set is 0-, C-, and then x-wave. For 
the lower frequency band, it is x-, 0-, and C-wave. There are five id- 
entical C-waves, plus five 0- and five x-waves for various combinations 
of Af and r g in each frequency band for one cycle of the program. 

We shall call these one-s waves FSK waves. 

Table 4.1 lists the phase changes between successive pulses in a 
carrier of all the FSK waves in the program. The 0-wave, x-wave, and 
C-wave conditions depend on the proper choice of t and A-f . 
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Table 4.1 


THE RELATIVE PHASE BETWEEN SUCCESSIVE PULSES AT A 
PARTICULAR CARRIER FREQUENCY FOR VARIOUS FSK WAVES 


Af (Hz) 
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There are other versions of the program based on the same principles 

with slight changes in the sequence as well as of the values of Af and 

X . In this report, those cases in which 2xAfr equals a multiple of 
s s 

rt radians will be discussed. Other versions of the program will not be 

presented here. Whenever examples are illustrated, the values of Af and 

x will be indicated, 
s 

A one-s FSK wave can be decomposed into its Fourier components which 
are discrete "spikes" in the frequency domain. The bandwidth of the 
spikes is about 1 Hz . In the time domain, each spike corresponds to a 
sine wave. Thus, a FSK wave can be viewed as the sum of many sine waves 
at various frequencies, phases, and amplitudes. All the sine waves last 
one second. The reason why the sidebands last as long as the FSK wave 
continues and the detailed derivation of the sidebands of FSK waves are 
discussed in Appendix C. 

In f-t space, an FSK wave could be represented either as segments 
of sine waves alternating between two frequencies at regular intervals 
for one second or as several one-second long sine waves with proper am- 
plitudes, frequencies, and phases. These two representations are sche- 
matically illustrated in Fig. 4.3. It is interesting to note that the 
same waveform can be represented by two different displays in f-t space. 
This is because "frequency" and "time" are not independent quantities. 

Tiie information contained in a "point" in the time domain is spread out 
on every point in the frequency domain. To obtain precise and complete 
information for a single point in the time domain, we must have infor- 
mation over the entire frequency domain and vice versa. 

But, in practical data processing, we are not willing to wait an 
infinitely long time in order to obtain precise information on frequency. 
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Thus, we chop off a set of data in the time domain and estimate the fre- 
quency components within the time period TV/ . There is definitely an 
uncertainty of the estimated frequencies. It has been shown [e.g., 
Blackman and Tulcey, 1958] that, roughly speaking, the product of the 
uncertainty in frequency Af and in time At must be greater than unity. 

Considering h ^ AfAt as a basic cell in f-t space where h > 1, 
we still have a choice of either Af or At in practical data processing 
devices. There are no "points" but "cells" in f-t space. The "shape" 
of the "cell" can be arbitrary so long as the "area" of the "cell" is 
conserved. Choosing a long At corresponds to choosing high frequency 
resolution (small Af), while a short At provides good time resolution 
but with a large uncertainty in frequency. 

In analyzing the FSK waves, we can resolve the pulses by making the 

time window (TW) , the length of data which are processed in a single 

scan, shorter than the individual pulse T g • On the other hand, we can 

resolve the sidebands of the FSK waves by setting TW > 4-r . A data pro- 

“■ s 

cessor, which can "print" the estimated information in f-t space, esti- 
mates the intensities of the frequency components of the signals within 
the time window and "prints" the information on a vertical line specify- 
ing an instant in the time domain. This "instant" is somewhat arbitrary 
within the period of TW. In many data processors, the location of the 
line is set to represent the instant at the "center of gravity" of a 
weighted time window. The uncertainty in the time domain can be speci- 
fied by the "second moment" of this window and definitely depends on the 
choice of TW. After the processor has finished the first scan, it slides 
its time window to a new segment for the next scan. The estimated infor- 
mation is "printed" on another line corresponding to another instant . 
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The process continues for the third scans, fourth scans, and so on. 
Usually, an appreciable portion of the data in successive scans are 
overlapped. It is important to point out that the sliding "distance" 
in the time scale between successive scans may be much smaller than the 
time window. But, a smaller sliding distance in the time scale does not 
imply a better time resolution. 

As illustrated in Fig. 4.3, the processor sees single frequency 
signals except at the edges of the segments when TWl is chosen. The 
signals are resolved as consisting of square waves alternating at two 
carrier frequencies. When TW2 is chosen, the processor resolves sev- 
eral spikes in frequency in each scan. Therefore, the resolved signals 
consist of several horizontal lines repre seating sinusoidal waves at 
various frequencies. Thus , depending on the length of the time window, 
the processor shows either segments of sine waves or continuous side- 
bands. 

An important question is raised : Do the energetic electrons that 
interact with the FSK wave in the magnetosphere see the wave as segments 
of sine waves or continuous sidebands? An answer is given in the next 
section. It has led us to develop a method of remotely estimating the 
length of the interaction region in the magnetosphere. 

C. Estimation of the Length of the Interaction Region in the 

Magnetosphere 

In the magnetosphere, the interaction between a WM wave and the 
energetic electrons resonant with it is thought to take place within a 
limited region, usually around the equator. When electrons move away 
from this region, the inhomogeneity in the geomagnetic field destroys 
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the resonance between the electrons and the WM wave, and there will be 
no average energy exchange between the WM wave and the electrons . Since 
there is no '’clear-cut" type of boundary separating the region in which 
interaction occurs from where there is no interaction, there must be an 
essential ai’bitrariness about every definition of the length of the in- 
teraction region (IR) . For example, Helliwell [1967] defined an inter- 
action region as the region over which the unperturbed phase angle be- 
tween the WM wave and a resonant electron remains within ± jt . Inan [1977] 
defines an interaction region as the distance between the equator to the 
point where the perturbed phase angle between the wave and an electron 
becomes rt-radian and shows that the length of the interaction region not 
only depends on the wave intensity, the electron pitch angle, but is also 
closely related to the initial phase angle between the electron and the 
wave . For a small wave amplitude , these two definitions differ by a fac- 
tor of two . These two definitions refer to the region within which elec- 
trons are scattered effectively. Whether the VLF wave that interacts 
with the electrons is amplified or not depends on the details of the 
electron distribution in the IR [Nunn, 1974; Dysthe, 1971]. An IR de- 
fined in these ways may then be categorized as a particle interaction 
region (PIR) , suggested by Helliwell [1978], 

It is believed that waves are amplified in the interaction region. 
There must be spatial variation of the wave intensity from one end of 
the IR to the other end in the steady state. Another definition of IR 
can then be stated as the region between the location where the wave 
intensity starts to increase spatially and to the location where the 
wave intensity becomes uniform spatially. The length of the IR defined 
in this way can be measured by satellites moving along a duct in the 
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magnetosphere. This IR is categorised as a wave interaction region (WIR) 
[Helliwell, 1978], 

The length of the IR we shall discuss in this section is character- 
ized by the averaged memory time of the resonant electrons. It must he 
related to Inch's definition of IR, that is. characterised by the time 
that the electrons remain trapped by the wave, The method is to find 

the critical pulse length, T eo> of the ESK waves . For the ESICwave 

T s ' GSS tIlan T sc’ thc electrons can resolve the sidebands of 
the waves . For a wave with ^ greater than the electrons can 

only resolve pulses of the FSIC wave. Therefore, ^ must be related 
to the electron's memory time. This effect provides a basis for measur- 
ihg the average electron's memory time in VLF wave-particle interactions 
in the magnetosphere. By assuming that electrons can remember only the 
wave that they have encountered during the interaction, the length of 
the IR has been estimated. 

Figure 4.4 shows the histogram obtained by counting the number of 
the FSIC waves which trigger emissions vs the values of Ts of the waves 

from the data between 1500 to 1525 TO on 11 October 1974. The program 

described in Section B was transmitted during this period. The abscissa 
indicates the values of Ts and the ordinate shows the number of the 
ISK waves which trigger emissions. Therefore, the vertical amis could 
be thought of as indicating wave triggering ability. The 0- and *-wave 
conditions are also indicated . The o-waves show up five times more fre- 
quently than any other waves in the transmitter program. The number of 

the c waves that trigger emissions has been divided by five in the his- 
togram . 
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The histogram shows that whether the IBS wave is in D- or *-=ondi- 

tion is immaterial to its triggering ability. The parameter on which 

the wave triggering ability depends is the pulse length V The waves 

with T s = WO msec, the longest pulse length among the five sets, show 

an emission triggering ability as high as that of C-waves and are much 

higher than the others. Among the other four groups of FSK waves with 

T = 10 , 20, 30, and 50 ms, the triggering ability decreases as T s M- 

s 

creases. 

This can be explained by two hypotheses: 

(1) Electrons can only resolve the segment features of the 
PSK waves when T s > 100 msec and the sideband features 

of the waves when < 50 msec. 

(2) There are wave-wave suppression effects among the side- 
bands . 

As indicated before, an FSK wave censists of two sets of EF pulse 
trains at carriers f, and f,. respectively. The electrons encoun- 
tering more than four pulses of an PSK wave during the interaction re- 
solve the sidebands. Each sideband can perturb and organise electrons 
within a finite V, range centered at the corresponding resonance ve- 
locity. men the frequency spacings between the sidebands are large, 

their perturbed ranges of electrons in V (| do not overlap . The Side- 
bands organise their own groups of electrons, and there is no mutual 
interaction between the sidebands. As the frequency spacings decrease, 
the perturbed ranges may overlap one another. The closer the spacings 
in frequency, the more the overlap. The chance of the electrons being 
organised by any one sideband coherently will be decreased, and hence 
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the wave triggering ability is reduced. This is a mutual suppression 
effect . 

rn Appendix C, we show that the frequency spacings between sidebands 

of a FSK wave are closely related to the value of Roughly speaking, 

the spacings are inversely proportional to t and independent of whether 

s 

the wave is in 0- or ^-conditions . As t s increases from 10 to 50 ms, 

the frequency spacings between sidebands are reduced from 50 to 10 Hz . 

There xs more mutual suppression among the sidebands as their frequency 

spacings become smaller. This explains why the triggering ability of the 

FSK waves decreases as t s increases from 10 to 50 ms, regardless of 

whether they are 0-waves or rt- waves . 

An explanation for the observation that FSK waves with T = 100 ms 

s 

can trigger emissions as often as the C-waves is postulated. The elec- 
trons have resolved the segment feature of the waves. Since t s is long, 
an electron may only encounter a single pulse and leave the interaction 
region before seeing a second pulse at the same carrier. Various elec- 
trons may see different portions of a pulse or different pulses. They 
will respond to individual pulses accordingly. By noting that a 100 ms 
pulse can trigger emissions by itself, we postulate that the emissions 
induced by the FSK waves with t s = 100 msec are triggered by individual 
pulses. A 100 ms pulse is long enough to produce a "self-sustaining” 
oscillation in the "wave -particle system" in the magnetosphere in this 
period. Therefore, it can trigger emissions as often as a longer pulse. 
The emissions triggered by short pulses may have different forms than 
those triggered by long pulses [e.g., Helliwell and Katsufrakis , 1974]. 

The former tends to be weaker, too. 
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In Fig. 4.5, three FSK waves with t =10, 30, and 100 ms are il- 
lustrated. These examples are taken from the data from which the histo- 
gram in Fig. 4.4 is constructed. The spectrograms in the upper panel are 
analyzed by setting the time window TW of the data analyzer to 100 ms . 

The analyzer can clearly resolve the sidebands of the FSK wave with T g = 

10 msec but not of those with t =30 and 100 msec . The spectrograms 

s 

in the lower panel are obtained by setting TW to 250 msec. Tae sidebands 

of the FSK waves with t = 10 and 30 msec but not those with t = 100 

s 53 

ms can be resolved. It is clear that emissions developed from the FSK 

wave with t = 100 ms are triggered by individual pulses . The emis- 
s 

sions developed from the FSK wave with 7 g = 10 ms are triggered by 
sidebands. The 4.6 kHz sideband triggers a rising emission at its end. 

The 4.75 kHz sideband triggers an emission that "steps up" to the 4.8kHz 
sideband at about the 400th msec and then forms a rising tone. Subse- 
quently, it is observed that the FSK wave with T g — 30 ms exhibits an 
interesting modulation pattern in the upper panel. A rising emission 
appears to develop from the upper cutoff frequency of the pattern. How- 
ever, the spectrogram with better frequency resolution shown in the lower 
panel indicates clearly that the modulation pattern is formed by the am- 
plified sidebands and that the rising emission is triggered by the side- 
band at 4.6 kHz . 

By Figs. 4,4 and 4.5, we have shown that electrons have a finite 

"memory time" during interactions with the FSK waves in the magnetosphere . 

The memory time is long enough for electrons to resolve the sidebands of 

the FSK waves with r <50 msec and short enough to respond to individ- 

s — 

ual pulses of the FSK waves with = 100 ms . 
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Fig. 4.6. ILLUSTRATION OF THE APPARENT INTERACTION TIME (AIT). 
The upper panel shows that an electron enters the interaction 
region, encountering a pulse traveling in the opposite direc- 
tion. The lower panel shows the situation where the .electron 
is leaving the interaction region. 



that must be met for the electrons to resolve the sidebands . Electrons 

can resolve sidebands of the FSK waves with t < 50 ms. Thus, 200 msec 

s — 1 

appeal’s to be a reasonable lower bound for AIT. 

Prom the electron resonance condition and WM dispersion relation, 
it can be shown [e.g., Helliwell, 1965] 


where 


V - 2v (1 - A) 
g P 

v „ “ T o = 

Vp = CA.C1 - A)] 1/2 


(4.4a) 


(4.4b) 


A = f/f H 

Therefore , 



(4.5) 


Using typical parameters at L « 4, 


f « 180 kHz 
P 

f Tr « 13 kHz 

f = 1/3 f 
u 

and assuming AIT = 200 msec, we find that 

L' > 2000 Km 


100 



Furthermore, we know that the resonant energetic electrons do not 

resolve the sidebands when t > 100 ms . This fact sets an upper bound 

s — 

on . We conclude that 

4000 Km > L r > 2000 Km ‘ 


D. Sideband Triggerings and Features of Coupling 

Multiple waves with small frequency spacings in between, propagating 
in the whistler mode (WM) in the magnetosphere , can intex’act with one 
another. A WM can interact with energetic electrons through cyclotron 
resonance, perturbing and organizing electrons contained within a small 
range of V centered at the cyclotron resonance velocity. Hie corre- 
sponding frequency range is termed the "coherence bandwidth." As the wave 
intensity increases , the perturbed V range grows and so does the co- 
herence bandwidth . Multiple waves often organize different groups of 
electrons at various ranges when the wave frequencies are largely 

separated. On the other hand, when the fx-equency spacings are so small 
that the pertui'bed x’anges of electrons overlap, electrons ox’ganized 

by one wave can also exchange enei-gy with another. Wave-wave intei-actions 
thus occur through these electi'ons . 

In this section, we shall show more examples of sideband triggering 
and then examples of enei’gy coupling between sidebands. 

Figure 4,7 shows an example of sideband triggering. The triggering 
wave is a rt-wave with r =10 ms and Af = 150 Hz . The corresponding 

'S' 

sidebands are spaced in frequency by 50 Hz. The calculated unamplified 
intensities of the sidebands are shown on the right. As noted in Appendix 
C, a FSK wave in a ^-condition is a carrier-suppressed signal. The two 
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Fig. 4.7. A TYPICAL EXAMPLE OF SIDEBAND TRIGGERINGS. Sidebands can be 
amplified and trigger emissions, just like any other constant fre- 
quency WM signal. 
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Fig. 4.8. MORE EXAMPLES OF SIDEBAND TRIGGERINGS. The panels showing the 
same data are obtained by successively reducing the analyzer's gain. 

The amplitude information is roughly obtained by this method. 
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becomes the only coherent signal near 4.02 kHz capable of organizing 
electrons efficiently. There are no other signals to suppress its growth. 
Therefore , the emission grows . 

The existence of an active PLH at 4.02 kHz is postulated because 
the emission frequency remains almost constant for about 500 msec. Fur- 
thermore, there are a few indications (not shown) of the existence of 
the 4.02 kHz PLH within half an hour, giving support to this interpreta- 
tion. 

This example indicates that the frequency spacings among the side- 
bands limit their growth . The coherence bandwidth in this particular 
case is greater than 50 Hz and depends on the intensity of the signal. 

It also illustrates two features of sideband coupling. The first 
feature in which energy steps up or down discretely in frequency is il- 
lustrated by the coupling between the 3.975 and the 4.025 kHz sidebands 
at about the 400th msec and that between the 4.075 and the 4.025 kHz 
sideband at about the 950th msec. The second feature in which energy is 
transferred between sidebands by emissions is illustrated by the couplings 
between the 4.025 and the 4.075 kHz sidebands at about the 600th msec. 
These two features of sideband couplings have been observed many times. 

E. Discussion 

The key to an explanation of energy coupling between two waves at 
different frequencies involves those electrons that can interact with 
both waves in the magnetosphere. An electron perturbed by a TO wave 
oscillates in V ( around the cyclotron resonance velocity V . The 
range of excursion in V (| depends on the wave intensity, the electron 
pitch angle, as well as the instantaneous phase angle between the electron 
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and the we when the electron is exactly at v = V q (see Appendix B). 
Various electrons in resonance with the wave may have different excursion 
range in V ( . As far as the wave is concerned, it can perturb and orga- 
nize electrons which are in a finite V range centered at the resonance 
velocity V o . We assume the perturbed V ( range to be the trapping range 
AV t* The corresponding frequency range is the coherence bandwidth. Two 
waves with a large frequency difference organize electrons in different 
V ll ran £ es • The waves are mutually independent. When the frequency dif- 
ference is small, the perturbed V (j ranges overlap. Electrons in the 
overlapped range can exchange energy with both waves . They serve as means 
of energy coupling between waves . 

According to the homogeneous model in which a constant wave is as- 
sumed (see Appendix B) , the perturbed V range AV ± is proportional to 
the square root of the we intensity. The wave intensity can be "mea- 
sured” by estimating the size of AV t . We have estimated the size of 

AV t by finding the critical frequency difference within which two waves 
interact „ 

It is important to point out that in reality this homogeneous model 
may not be adequate because 

(1) the inhomogeneity of the geomagnetic field changes the 

size of AV, , and 
t 

(2) the wave intensity in the interaction region is not 
uniform . 

To the author's knowledge , there is no published relationship be- 
tween the size of AV t and wave intensity (even assuming a constant 
wave) in inhomogeneous models. 
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According to Dysthe 1 s simple model C1D71] in which a constant ex— 
terna! torque is assumed to account for the inhomogeneity in an interac- 
tion region (IR) , the size of AV t is reduced as the IR moves away from 
the equator. Thus, to produce the same AV^ , a stronger wave is required 
in the inhomogeneous model than in the homogeneous model . Since WPI is 
believed to occur around the equatorial region where the inhomogeneity 
is small, we shall use the homogeneous model to estimate the wave inten- 
sity in the IR. The wave intensity in an inhomogeneous IR must be greater 
than the estimated value. Moreover, the actual wave intensity is not 
uniform in the IR. Therefore, the estimated value can only account for 
the order of magnitude of the average intensity in the IR. 

It is interesting to visualize the coherence bandwidth from the point 
view of electrons . A constant frequency WM wave interacts with energetic 

electrons in a finite IR, causing the electrons to oscillate in V . The 

II 

Doppler-shifted frequency of the wave as seen by an electron with a small 

velocity excursion AV ( is not constant. The wave has a finite bandwidth 

kAV (j as far as the electron is concerned. This bandwidth is a "resonance 

bandwidth produced by WPI. Electrons with various AV see different 

bandwidths (kAV ) associated with the wave. The electrons with AV = 

II 

AV t see the lar eest bandwidth which we call the coherence bandwidth. The 
coherence bandwidth is in the order of 50 Hz or more for most of the co- 
herent m waves in the magnetosphere. 

When an electron encounters two waves in the magnetosphere, its mo- 
tion becomes more complicated. Suppose that the two waves at f and f 

1 2 * 

traveling with the phase velocities V and V , have intensities 

P*** P" 

^wl’V 5 and (E w2 ,B w2 ) » respectively, where >± 2 , V ± >V £ , and 
B wl > B w2* total wave magnetic field B^, as shown in Fig. 4.10, 
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Fig. 4.10. AN ELECTRON MOTION IN TWO WAVES. 
The frame is moving with the phase velocity 
of wave 1 . The magnetic field of wave 2 
rotates with an angular frequency that 
equals the frequency difference of these 
two waves. 
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is a vector sum of B „ and B „ and varies with time in the frame 
moving with V pl * We have assumed that l V pi“ V p 2 l is mucil smaller 
than V pl and V 2 so that the electric field of wave 2 is negligible. 
In that frame, disappears, B ^ is static [e .g . , Matsumoto , 1972] , 

and B „ rotates with a frequency 0, where 9 equals -f |. 1/0 


corresponds to a time constant that equals the rotation period of B 
and is called i.e., 


w2 


T I = 1/0 


(4.6) 


This time constant t shall be the same in the lab frame. The 
reason we derive this quantity in a wave frame is to "fix" one wave in a 
frame so that the relative rotation between the field components of these 
two waves can be visualized easily. 

The detailed calculation of energy coupling between two waves through 
the resonance electrons is still being investigated. Meanwhile, we can 
still argue qualitatively for a rough criteria that determines whether 
the waves will mutually interact or not. It is believed [e.g., Helliwell, 
1967; Helliwell and Crystal, 1973] that the bunching time is a char- 

acteristic time of energy exchange between a constant amplitude wave and 
its resonant electrons. When f ~f is large, the electrons in reso- 
nance with wave 1 can not effectively interact with wave 2. In one bunch- 
ing period , has rotated more than one cycle . its effect on an elec- 

tron which is in resonance with wave 1 has been smeared out. Similarly, 
electrons in resonance with wave 2 cannot effectively interact with wave 
1. Thus, these two waves are organizing different sets of electrons and 
do not mutually interact. One could specify a criterion under which two 
waves are mutually independent as 
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T b > Tj C4.7) 

Accoi’ding to Eqs. (B.16) and (4.6), we can rewrite the criterion as 

AfcB < l f l ‘ f 2 i ' (4.8) 

It states that, when the frequency difference is greater than the coher- 
ence bandwidth, two waves do not interact. Thus, by finding the critical 
frequency difference within which two waves interact, we have found the 
coherence bandwidth The wave intensity can then be estimated by 

Eq. (B.16) . 

Equation (4,8) is a rough criterion. We have assiqged to be 

small enough so that it does not alter the motions of electrons which are 
in resonance with wave 1. This 4s a "small perturbation" assumption. 

When B w2 is not small, a criterion that determines whether two waves 
interact or not must involve both B 1 and B _ , One such criterion is 


2 tAC CBl * < l f l " f 2 


(4.9a) 


where Af CBl and A£ cb2 are the corresponding coherence bandwidths of 
wave 1 and wave 2. It implies, as shown in Fig. 4.11, that the perturbed 
V ranges associated with the two waves in the electron distribution 
overlap unless * 

| 0*,,, 4 AV ti ,) < |V 01 - V Q2 J (4.9b) 


where and Vq 2 are the resonance velocities associated with wave 


1 and wave 2, respectively. 
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DIFFERENCE. The purturbed V„ ranges are not overlapped. 
There she!! be no interaction between these two waves. 


From Eq. (4.9a), it can bo shorn that there is a factor of 4 in the 
estimated intensity of wave 1, B^, between the cases in which = 

B w2 and the ones in which B wl » B^. Since wo are only interested in 
the order of magnitudes of the wave intensities in the IB, Eq. (4.S) ap- 
pears to be suitable for our purpose. We shall bear in mind, however, 
that the estimated value could be off by a factor as largo as 4. 

The excursion range AV (] for various electrons in resonance with a 
wave are different. It seems reasonable to average these excursion ranges 
It can be shown by Eq. (B.4b) that 


AV = 2 


n v 

2 ■■ (1 - cos ty_) 


Iv 


R 


1/2 


(4.10) 


whore is the instantaneous phase angle between the wave and an elec- 

exactly equals V . Using the approximation for trapped 


R 

tron when V 


electrons 


V ~ V tan a 
i o 


(4.11) 


we have 


■Of 

4/2 | . R 
IcAV = 4 (kfl V tan a) |sm ~ 
It w o 


( 4 . 12 ) 


Assuming that electrons initially are uniformly distributed in the phase 
angle, we can show that AV u averaged over is 


IcAV 


* 


S 1/2 1/2 

- (kfl V ) ' tan a 
jt wo 


(4.13) 


R 


Thus , we have an 


averaged coherence bandwidth (over the phase angle) 
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AF 


CB 




R 


4 <“wV 1/2 

it 


tan 


1/2 


a 


(4.14) 


Equation (4.14) reveals that Af CB | is dependent on the electron 
pitch angles. In the magnetosphere, waves can interact with electrons 
at various patch angles. The electron distribution as a function of en- 
ergy and pitch angle in the magnetosphere can be chosen as [e.g., Lietn- 
ohn, 1967] 

f (W,a) « W“ P sin m a (A i S \ 


This model has several advantages. It is mathematically simple. 
The energy variation is reasonable and can be fitted into the particle 
data [S chi eld and Prank, 1970]. 

Equation (4.15) can be written as 


f( VV a) « 


• ni in 

sin a sin a. 

— OC - 


2 2 
V + V 


,_2p 2p 

V sec a 
o 


(4.16) 


i . e . , 


f(0i) « sin m ex cqs^ cs 


2p 


(4.17) 


We have used the resonance condition V (j « V q and V. ~ V q tan a. Fig- 
ure 4.12 shows the distribution fraction f(a) for various values of m 
and p. It is noticed that the peak tends to move to high pitch angles 
as m increases and tends to move to low pitch angles as p increases. 
The condition m = 0 implies an isotropic pitch angle distribution at a 
given energy level. It does not imply uniform distribution at a given 
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Fig. 4.12. PLOTS OP •pOy'l — .,4 2p 

VARIOUS VALUES OF m AND p. The T “ a t wh’ 0F “ TOH 

value is at q ; = tan" 1 m/2p . max 1 h h f Ilas a maxin 



We shall average the coherence bandwidth again over a pitch angle 
distribution. From Eq. (4.14), 


where 


08 = “CB 


* 


IS M/ b r(«,p) 


r(m,p) « f 0 ^ tan 1//2 a g(o0 da 


(4.18) 


grCce) = 


TTl 2p 

sin d cos a 


f 


^ Sin m a cos 2p a da 


A loss cone 


“l = 


- V 

A -1 o 
% = cos “ 


Figure 4.13 shows the values of T as functions of m and p for a 
case in which a L = 7° and 0^ = S5°. It is on the order of unity. 

Equation (4.18) is the relationship between the coherence bandwidth 
and the wave intensity in the magnetosphere. We have a qualitative cri- 
teria for deciding whether two waves with a frequency difference Af can 
mutually interact or not. men Af > CB, there is no mutual interac- 
tion. When Af < CB, they will interact. Thus, by estimating the crit- 
ical frequency difference Af c within which two waves interact and using 


Eq. (4.18) and the condition 


Af * CB 
c 


we have been able to measure the average wave intensity in IR. For sim- 
plicity, we shall refer to CB as "coherence bandwidth" hereafter. 
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Fig. 4.13. PLOTS OF r(m,p) FOR THE CASE WHERE 
OL - 7 AND a ~ S7°. 
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•F • LISP Program 

From Section 0, we have found that sidebands with frequency differ- 
enoe less than 50 He can mutually interact. We have also found in many 
occasions that sidebands with a 50 Hs frequency difference can grow and 

trigger emissions independently. It appears that 50 He could be a good 
measure of coherence bandwidth. 

mere are several transmitter programs designed specifically to 
measure the coherence bandwidth. One such program is the Line Spectrum 
(LISP) Program. The LISP program, she™ in Fig. 4.14, consists of ESK 

waves with 1,-10 msec. The frequency difference between the pulse 

train at the upper frequency and the one at lower frequency is switched 
between 100 and 200 Hs every minute. Ihe upper frequency remains 

the same during the transitions . As a result of this modulation, the 

spectrum of the wave is switched between the two shorn in Fig. 4.15. The 

sidebands are spaced at multiples of 50 T+ . .. 

xp-i.es 01 ou Hz. it is noticed that the side- 

band at f HI is always at the same amplitude and continuous in phase, 

while its nearby sidebands change either amplitudes or phases during the 
transition. 

If the coherence bandwidth of the sideband at f^ is i ess than 50 

Ha, then its behavior is independent of the changes in nearby sidebands. 

On the other hand, if its coherence bandwidth is greater than 50 Hs, it 

interacts with other sidebands at nearby frequencies. Therefore, alter- 

mg the nearby sidebands should affect its » ow +i, „ . 

u di-iecr its growth. By examining the am- 

plitude of the sideband at f durine- f ra u Sn -+i n 

HI ing transit tons , we can determine 

whether the CB is greater than 50 Hz or not. 

Figure 4.15 shews an example „f this program. The upper panel shows 
the spectrogram in which the time and the frequency resolutions are about 
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tt/2 

-tt/2 


I ■ I 


'LO 


‘HI 


Fig. 4.15. THE SPECTRUM OF AH FSK WAVE WITH 
ms AND THAT WITH Af = 200 Hz AND 
both spectra have the same amplitude and phase, 
spectrum. 


T s = 10 ms . 


Af = 100 Hz AND tt s » 10 


The spike at f HI in 
0 is the phase of the 
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Fig. 4.16. AN FXAMPLE OF THE LISP PROGRAM. The dynamic spectrum has high frequency resolution 
(~ 2 Hz). The stimulated sidebands have been resolved,. 


350 msec and 3 Hz, respectively. The middle and lower panels depict the 
amplitude variations of the 4.35 and the 4.15 kHz sidebands. The band- 
width of the filter through which the amplitude charts are made is about 
50 Hz . It is important to point out that the records have a compressed 
time scale. The regular fluctuations appearing on the signals in the 
spectrogram are due to asynchronization between the recording and play- 
back recorders . 

There are Pill's on the spectrogram. The strong signal around 4.26 
kHz drifting slowly in frequency at the rate of a few hertz per minute 
is a local induction line (71st PLH). The 4.25 kHz sideband appearing 
every other minute is observable, with difficulty, only in the first, 

the third, and the seventh minute, where the induction line is above 
4»255 kHz • ...... . 

There are "stimulated sidebands" about 7 Hz away from the parent 
sidebands. A good example can be observed in the eighth minute. The 
expanded dynamic spectrum of the eighth minute is shown on Fig. 4.17. 

The 4.20 kHz sideband has stimulated two sidebands spaced at about 7 Hz 
above and below the parent sideband. The 4.15 and 4.35 kHz sidebands 
also stimulate sidebands in this period. It is observed that the stim- 
ulated sidebands occur only when the parent signals are strong. It has 
been reported previously that a constant frequency signal can trigger 
sidebands [Chang and Park, 1977]. Stimulated sidebands are believed to 
be a phenomenon caused by nonlinear wave-particle interactions in the 
magnetosphere . 

The 4.35 kHz sideband is at f^, the frequency at which both the 
sideband amplitude and the phase remain unchanged at the transitions. 

The slow changes in amplitude shown on the middle panel are due to the 
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IN THE EIG^T^TE 1 IN^ G ? F 4 !i6. DYNAAIIC SPECTR ™ 
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SlOB te " POral *>■ «■* magnetospherie conditions. The transi- 

tions occur at ~2.3 seconds after the minute marks. From such compressed 
time-scale records, it is hard to locate precisely the time of transi- 
tions. However, the lower pane!, showing the amplitude of the 4.15 kHs 
Sidehand that is "switched on" every other minute, provides definite in- 
dications of the time of the transitions, it is noticed that, at the 
transitions at the beginnings of the third, fourth, and tenth minute, 
there are sudden amplitude variations on the 4.35 kHa sidehand. The 
variations can he as large as 3 dB. from the studies of the expanded 
records, it is clear that there are no sudden amplitude variations dur- 
ing the rest of the transitions. The expanded records during the tran- 
sitions at the beginnings of the second, fourth, eighth, and tenth min- 
ntes are shown in Figs. 4.13, 4.13, 4.20, and 4.21, respectively. Ko te 
that the 0-dB levels in various amplitude records do not refer to the 
same absolute amplitude. The arrows indicate the time of transitions. 

The frequency resolution of these spectrograms is different from the one 

shown on Fig. 4.16 and is about 7 Hz m, p -in , 

H . The sidebands stimulated by the 

4.35 kHz sideband cannot be resolved. 

it is observed from Fig. 4.16 that only when the 4.35 kHa sidehand 
is at high intensity can the transitions affect its amplitude. This sug- 
gests that the intensity of the 4.35 kHa sidehand in the X* is such that 
its coherence bandwidth is about 50 Ha and varies slowly with time. Hs- 
ihg the available whistler data, we have found that the path is at h„4 

^ %“ 40 °™' 3 - Assuming that p = 4 and m = 1 for the electron 
distribution in the magnetosphere and a CB equal to SO Ha, we have cal- 
eulated the wave intensity in the IB by Sq. (4.13) to be on the order of 
2.5 »7. Ffom the argument given in the previous seetion, there could he 
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Fig. 4.18. EXPANDED TIME SCALE 
SECOND MINUTE IN FIG. 4.16. I 
same as shown in Fig. 4.16. 
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Fig. 4.21. EXPANDED TIME SCALE OF THE TRANSITION AT 
THE TENTH MINUTE IN FIG. 4.16. O-dB levels are 
not the same as shown in Fig. 4.16. 
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an as ^ as a factor of 4 in t* estimated wave intensity. T hus 

the estimated wave intensity is somewhere between 3.5 to 10 my, to gooa ' 

agreement with satellite measurements [Heyborno, ,1986, Burtis , 1974 • I„an 

et ax, 1977J . estimation is based on the assumption that the’ in- 

teraction occurs at the eauntm- M ,, . 

equator . If the interaction occurs off equator 

where the inhomogeneity of the geomagnetic field h. 

g 1C field becomes important, the 

estimated wave intensity would be larger. 

It is interesting to note that the intensity of the 4.35 ** 

band increases at the transitions in the third and • „ 

thud and ninth minutes where 

" 4 ‘ 35 ^ dem ’ eaSM « in the fourth min- 

ute where the 4.15 Ms sideband is turned on Tt • 

_ • It is suggested by Kelli- 

X Cprxvate communication, 197SJ that the 4.4 kHz sideband « t 60 

He in Pig. 4.15, may coupie energy to the 4.35 * sideband at th! tran- 

sit ions where the intensity of the 4 35 h-Wt - . 

the 4.35 kHz sideband is high. Bie dyna- 

spectrum in Fig. 4.1s indicates that the 4.30 kHz C f - 50 He) 

sideband has not been amplified significantly in « HI 

s mcantlj- in the magnetosphere. At 

the transitions where the 4.15 kHz 1 , 

sideband is turned off, the 

input intensity of the 4.4 kHz sideband is reduced b, 10 

i sauced by 10 percent, it 

should therefore be amnl •? i 

... eSS ' and the corres Pond ing coherence band- 

width should decrease 

decrease. Therefore, the intensity of the 4.35 kHz sideband 
increases, consistent with the observation. 

G ‘ Summary and Conclusion 

in this chapter , we have shew, that the PSK waves can be viewed as 
several constant frequency signals, in the magnetosphere, electrons can 
resolve these sidebands *en the modulation period in the signal is small v 
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compared to the interaction time between the wave and the electrons. 
Otherwise, the electrons can only resolve pulses of the FSK waves. This 
xs the basis of the method we used in estimating the length of the in- 
teraction region CXR) . It has been found that the IR lies between 2000 
and 4000 Ion. 

Transmitting FSK waves is equivalent to injecting multiple waves 
into the magnetosphere. The frequency difference is controllable by 
varying the pulse length. We have found that waves with frequency spac- 
ing less than 50 Hz often show mutual interactions. Suppression and 
coupling among sidebands are often observed. We have identified two 
different features of energy couplings. In the first feature, energy 
of a sideband transfers to another sideband by a discrete jump in fre- 
quency. in the second feature, the energy coupling is accomplished by 
"small’' emissions that link one sideband to another. 

The interaction between two forced mode signals is attributed to 
the electrons that can exchange energy with both waves. Assuming that 
the trapping width is the width of V| , within which electrons can be 
organized effectively by a wave,.* have established the relationship 
of wave intensity and the critical frequency spacing within which two 
waves can interact. This critical frequency spacing is called the co- 
herence bandwidth and is determined by averaging over both the phase 
angle as well as the pitch angle, it is found that the coherence band- 
width is on the order of 50 Hz for the sidebands. The corresponding 
wave intensity is somewhere between 2.5 to 10 m 7 in a homogeneous model 
and could be larger in inhomogeneous models. 
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Chapter V 


SUMMARY AND RECOMMENDATIONS 

A. Introduction 

The purpose of this report has been to present an experimental in- 
vestigation of VLF wave-wave interaction in the magnetosphere. In this 
chapter, we shall summarize the contributions of this report and make 
suggestions for further investigation. The contributions in the study 
of pulse distortion, due to propagation in the magnetosphere, in the 
area of interaction between gap-induced emission and post-gap signal 
and in the field of sideband mutual interaction will be covered. The 
suggestions for further work will include the same three areas. 

B . Summary of the Present Work 

1 - Pulse Propagation in the Magnetosphere 

We have been performing experiments to investigate VLP wave- 
particle and wave-wave interactions in the magnetosphere by injecting 
VLF pulses from the Siple transmitter. These pulses are distorted be- 
fore arriving at the interaction region which is believed to be near the 
equator. The distortion caused by propagation through the dispersive 
medium may have stretched the pulses enough^to alter the details of WPI 
in the magnetosphere. It is therefore essential to know how the pulse 
is distorted due to propagation in the magnetosphere. 

As far as we know, the distortion of a VLF pulse propagating 
in the magnetosphere has not been calculated previously. 
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In Chapter II, we have described the method of analysis. 
Study of pulse propagation in the magnetosphere has provided the follow- 
ing results : 

(1) The group velocity calculated at the carrier frequency 
can accurately predict the arrival time of a pulse only 
when the time delay refers to the "centers/' not to the 
leading edges of the pulses . 

(2) The frequency of a distorted pulse is not constant over 
the duration of the pulse. The largest deviations occur 
at the front and/or the rear ends of the pulse. 

(3) The frequency at which a pulse lias minimum distortion at 
the equator is the "nose" frequency. 

(4) A discontinuity in the phase of a fixed amplitude sine 
wave introduced on the ground can develop into an ampli- 
tude gap when the wave arrives at the equator. The gap 
can be as large as 20 msec as a wave at 3 kHz propagating 
through a duct at L = 4. 

(5) A digital equalizer has been designed to compensate for 
the phase distortion. Computer simulations show that the 
equalizer can completely cancel the phase distortion. 

This device could be used either as a pre— or post— pro- 
cessor to remove the distortion developed in the magne- 
tosphere . 


2 . Gap-Induced Emissions 

The phenomenon of gap-induced emissions has been reported for 
the first time. A 10 msec gap in a tz’iggering wave can induce emissions 
which may then interact with the post-gap signal, causing suppression and 
entiainment effects. The interaction involves a natural mode and a 
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forced mode signal, and this is termed forced mode/natural mode wave- 
wave interaction (ENWWI) , 

Among the results of this study are the following: 

(1) This phenomenon suggests that the tail-end emissions 
start to develop within a time period less than 10 ms . 

(2) Experimental data show that, as the gap is opened up 
in time, we can "scan" the development of emissions. 

One of the results confirms an earlier observation 
that emissions always start with a small rise in fre- 
quency. 

(3) A 70 msec gap is long enough to permit a falling tone 
to fully develop . 

(4) A new feature of FNWWt has been discovered. The post- 
gap signal can "capture" a falling emission that is 90 
Hz below the transmitter signal, causing the emission 
to overshoot toward the positive frequency direction. 

The frequency difference between these two signals ap- 
pears to be an important factor in determining whether 
or not the interaction can occur . 


3 . Sideband Triggering 

Sideband triggering has been reported for the first time. It 
was an unexpected result from the studies of a transmitter program de- 
signed for other purposes. This discovery enabled us to study sideband 
interp.ctions . Suppressions and couplings between sidebands are commonly 
observed . We call this interaction SIDEBAND MUTUAL INTERACTION (SMI). 

From modulation theories, we know that discrete sidebands are 
created by periodic modulations on the signals. However, in the 


135 


magnetosphere, the energetic electrons respond only to certain sidebands 
The electrons can interact with the sidebands of the waves with small 
modulation periods but not with those generated by long modulation pe- 
riods . We have shown in Chapter iv that this observation enables us to 
estimate the length of the interaction region. 

Study of sideband triggerings has provided the following re- 
sults . 

CD The length of the interaction region has been estimated 
to lie between 2000 and 4000 km. 

(2) Sidebands, like any other constant frequency signals, 
can trigger emissions. 

(3) 50 Hz appears to be the critical frequency spacing within 
which sidebands mutually interact. 

(4) The coupling is attributed to the electrons that can 
interact with both waves. It is believed that a wave 
can only organize electrons within a finite range of 
V |( centered at the resonance velocity. Thus, the 
interaction must come from the electrons in the over- 
lapped 'V. range . 

(5) The LISP program suggests that 50 Hz is the order of 
the coherence bandwidth of sidebands. The sideband 
intensity in the interaction region has been estimated 
to be on the order of 2.5 ~ 10 my according to a homo- 
geneous model, in reasonable agreement with satellite 
measurements Cffeyborne, 1966; Burtis , 1974; Inan etal, 

1977].. 
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C. Suggestions for Future Work 


1 • Pulse Propa gation in the Magnetosphere 

The work can be extended to nonducted signals by combining a 
ray tracing program and the algorithm discussed in Chapter II. it is 
generally true that the dispersions for nonducted signals are greater 
than for ducted signals. Therefore, the propagation distortion shall 
be even more important in the nonducted case. 

It is necessary to use a pre-processor, or equalizer, on the 
transmitter end to compensate for the phase distortion developed from 
the path between the transmitter and the interaction region. Realiza- 
tion of this device appears to be feasible. 

2 . wvm 

It seems feasible to design a transmitter program to study 
the capture" feature of the FNWWI in a somewhat controllable way. The 
slopes of the frequency of the gap-induced emissions can not be con- 
trolled easily. But the frequency and phase of the post-gap signal can 
be changed. Thus, studies of this intei'action for various conditions 
on the frequency and phase differences between the falling emissions 
and the post-gap signals can be carried out. 

3. SMI 

More VLF wave injection experiments are needed to establish a 
relationship between the input wave intensity and the size of the 
coherence bandwidth. At present, we only have qualitative evidence , 
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suggesting that coherence bandwidth increases as the wave intensity be- 
comes stronger , 

It will be very useful to inject only two WE waves into the 
magnetosphere . The frequency difference will be varied from a few hertz 
to 200 Hz for various wave intensities . This experiment will give us 
quantitative information on the coherence bandwidth versus the input 
wave intensity. It may also tell us when the mutual suppressions occur 
and when the mutual couplings become predominant. A refined model of 
wave-wave interaction in the magnetosphere can then be constructed from 
this information. 

In the theoretical area, it is necessary to develop the studies 
of two-wave interactions , either by analytical methods or by computer 
simulations . Understanding the electron motions in two WM waves in a 
homogeneous model appears to be the first step in the solution to this 
problem. 
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Appendix A 


A COMPUTER CODE SIMULATING VLF PROPAGATION IN THE MAGNETOSPHERE 


This appendix lists one of the computer programs that simulate VLF 


pulse propagations in a duct in the magnetosphere. 


™ E p?pp° RTED PULSED IS TRANSLATED TO 500 HZ FOR DISPLAY 
C THE CARRjErfSEQSL^cJ^iS^H^^JARJlSo^iSE OF^Je JSlse' 

c 25 5?!!c!S!^5 TD * -■» 

c THERE ARE 6 DISPLAY OF THE PULSE SHAPE ALONG THE PATH FOR EACH FREQUENCY. 


II 


60 


110 


111 


61 


112 


COMMON B(1026) .BIC 10261 , CR ( 51 2 ) , Cl ( 512 ) , PD ( 51 2 ) 

REAL BI,CR,CI,PD,SI(126),CT (1025) 

REAL L»FQ,DPH,DDFP,NEQ 

INTEGER NF,SA,DE,NP,NL, INC65) ,NFR 

INTEGER N3 

EQUIVALENCE ( BI ( 1 ) , CT ( 1 ) ) 

CALL TPLTR 
CALL PIQCS( 20 ) 

WRITE (1.11) 

FORMAT ( * NO . OF FILES TO BE SKIPPED IN TAPE 82») 

READ (l.fc) NF 
CALL SKIPFCS2.NF) 

READ (5 , s)L,DE,NEQ,NP, F0,NT 
N=512 

WRIJE (6,60) L.DE.NEQ 

NI = 2^NT+1 X ' ,L " T ' F5,1 ’ 2X ' ,DE=, ' n ' 2x »'NEQ =’ *F10.2» ’PER CO 
NL=2« (NT+NP) 

READ ( 5 , s ) NFR.DF 
DO 110 1=1,1026 
B ( I ) = 0 

DO 111 I=NI , NL , A 
B(I) = 0 
B(I+1)=1 
B C 1 + 2 ) = 0 
B ( 1+3 ) =-l 

CALL CTA (10, B, IN, SI, 0) 

CALL CTA ( 1 0 , B , I N , S 1 , 3 ) 

DO 1011 K=1,NFR 
INT =17 

WRITE (6,61) F 0 , NP 

FORMAT (/IX, 'THE CARRIER FR EQUENCY= ' FI 0 . 2 , ' HZ ' , 2X . 

1 'THE PULSE LENGTH =',I3,' MS') ' 

DO 112 1=1,1026 
MM= 1 

BICI)=B(I) 

CALL CTA (10, BI, IN, SI, -3) 

BI(1025)=BI(1) 

N3 = l 

CALL PDIPO(L,FO,NEQ,DE»NP) 

CALL PLOPS (CT.N3) 

DO 1010 KI=1 ,5 
N3=N3+1 

CALL PDELY( L.DE.NEQ, INT, DPH , N, FO ,MM) 
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cz 


1010 

1011 


c 

c 


DDFP= ( INT“1 JkDPH* 180/3 .14159 

WRITE (6.62> DDFP TPAVELED FROM THE BOTTOM OF THE DUCT UP’ 

/?T?5E X ^rniSl! L AWA5iRmi A ^E E EQUATOR BY '.F10.*.* DEGREE* 1 

CAUL CTA C10iBI.IN.SI.-3) 

BIC1025)=BIC1) 

CALL PLOPS CCT.N3) 

int=int-a 

fo=fo+df 

CALL TPLTX 

STOP 

END 


SUBROUTINE PDELYCL.DE.NEQ.INT.DPH.N »FO >MM3 


THI5 SUBROUTIHE CREATES THE PHASE BEEAV BP EACH PREBUENC COMPONENTS 
FOR A PULSE PROPAGATING IN A DUCT. 


COMMON BC1G26).BIC1026)>CRC51H)»CIC512)»PDC512) 

SKI; FHC21>,SF,CF,NEQ,0M,L,LM,KE.HH,LM1 

INTEGER 1NT. IST2 . IST3 . N , Ml . J R » J I » BE 
Rl=6500 

RO=6AOO , 

Gl=9 • 8E-03 *(R0/R1>**2 

l THE BOTTOH OP THE MAGNETOSPHERE UP 

c TO THE EQUATER 
REQ=R0SL 

FPEQ=9E03 *5QRTCNEQ) 

RT= SQRT C(R0*L-R1)/R1> 

FAI 1= ATANCRT) 

SF1= SINCFAI1) 

CF1= COS t FAI 1 1 
IF CDE .EQ.l) GO TO 100 
IF C DE .EQ.2) GO TO 101 
IF CDE .EQ.3) GO TO 102 
IF CDE .EQ.A) GO TO 103 

100 CH=0.08 
CO=0.9 
CHE=0 . 02 
T = 1 600 

GO TO 200 

101 CH=0 . 08 
C0 = 0 . 9 
CHE=0 . 02 
T=320B 

GO TO 200 

102 CH = 0 . A 
C0 = 0 . 5 
CHE=0 . 1 
T = 1 6 0 0 

GO TO 200 

103 CH=0. A 
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200 


‘if- 

'T 


1000 

.!&>■ 


60 

61 

608 

69 

700 


C0 = 0.5 
CHE=0.1 
T=B00 
FAI = 0 



KE=1 . 38E-23 

HH=1 . 673 E-27 

HH= KE*T/CMHkG13#1.E-06 

HD=HH/16 

HHE=HH/6 . ^ _ 

ZEQ=Rl-< R1*S23/REQ-(0M**2/C 2*G1 3 3*( REQK*2- C R1*CF1 3 *H2 3 

EHEQ= EXPC-ZEQ/HH3 

EOEQ=EXP(-ZEQ/H03 

EHEEQ= EXP ( -ZEQ/HHE 3 

vde=ch*eheq+coseoeq+che*eheeq 


LH1» SQRTC 1/VDE ) 

TEMPI- SQRTC1+3SSF1HSF1) 

PI ( 2 1 3 = R0*L*CF1*TEMP1 
FPC21 3 = S<3RTCLN13*FPEQ 

FHC21 3= B.736E05*( 1/L/CF1/CF1 3**3*TEMP1 

DPH-FAI 1/20 

DO 1000 1=1.20 

SF = SINCFAI3 

CF = COS ( F AI 3 

R=RO*L«CF*CF 

Z=Rl~<Rl s *2)/R-tOMs*2/C2KG113*((R£CF3**2-(Rl#CF13**23 
EH=EXP (— Z/HH 3 
E0= EXPC-Z/HO) 

ERE=EXPC-Z/HHE) 

VNU=CH*EH+CDsEO+CHE*EHE 

LM= SORT ( VNU/VDE 3 

TEMP = SQRTC1+3«SF*SF3 

PI ( I 3 = R0SL*CF*TEMP 

FP ( 1 3 = SORT ( LM3SFPEQ 

FH< 1 3 = 8.736E05H(1/L/CF/CF3K*3*TEMP 

FAI = FAI+DPH 

IF CMM .GT. 1) GO TO 608 
DFAIl=FAIl«18Q/3. 16159 
DDPH=DPH si 80/3. 16159 
PRINT 60.DFAI1.DDPH 

FORMAT (1X,'DFAI1=',F7.2, ’DEGREE* . ,.ZX» ' I MORE . = ' . F7 . 2 , ' DEGREE ' 3 
WRITE <6,613 C FPC I 3 , FH ( 1 3 » PI < 1 3 » I = K5'Sl 3 
FORMAT < 3F12 .6 3 
CONTINUV 

IF CINT .LE. 173 GO TO 700 
WRITE (6,693 

FORMAT ( 1 INT IS TOO LARGE. INT SHOULD BE 1,3,5, UPTO 17'3 
STOP 

IST2=INT+1 

IST3=INT+2 

RIC2= FPCINT3*«2/FQ/<FHCINT3-F03 


DO 900 J=1,N 
FW=1000 

FI (J3=F0-FW/2+FW«J/N 
F=FICJ3 

RI2= FPCINT3*tf2/F/(FH<INT3-F3 

RRI = SQRTCRI2/RIC23 

FR I = SQRT ( RR 1 3 

DO 890 1= INT, 21 

GT = F# C FH C I J — F J 

G=FP(I3/SQRT(GT3«P1(I)/3E05 
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o o n o o on 


S9D 

891 

898 


695 


900 


TY(I)=DPH#G/3 
TT=TY(INT)+TYC21 J 
do 891 i=isra,2o,a 
TT=TT+A*TY( I ) 

DO 892 I=IST3,19,2 
TT=TT+2«TY( I ) 

TYL=TT 
CYL=F*TYL 
HI = CYL 
DIF=CYL-M1 

IF (DIF .LE. 0.5) GO TO 895 
DIF=DIF-1 

FASE=DIF«2«3 . I<tl 59 
PD ( J ) = 36 OHDI F 
SN=5IN(FASE) 

CS=C05 { FASE ) 

JR=2XJ+1 
JI=2*J+2 
ARE=BC JR) 

AIM=B(JI) 

BI< JI )=£-AREkSN+AIM*CS)/TRI 
BI(JR)=(ARE*CS+AIN*5N )/FRI 

CR(J)= BKJR) 

CICJ)=BICJ1) 

BXC1)=B(1) 

B1 ( 2 ) = BC2) 

RETURN 

END 


SUBROUTINE PLDPS(CT,NCL) 


THIS SUBROUTINE PLOTS SIX VLF PULSES AS FUNCTIONS OF TINE 


690 

691 


700 

701 


REAL CTC1025) 

INTEGER NCI. 

IFCNCL .GT. 1) GO TO 690 

CALL RAML 

XO=-250 

Y0=2 

GO TO 691 
X0 = Q 


Y0 = 3 

CONTINUE 

CALL SCALE(0.01,0.5,X0, YO) 
CALL GRID ( 1 . 0. , 0 . , 1 . f 1 ) 
CALL GRID (3,0. ,0. ,1.,1) 
CALL GRID (1,-2 0. ,0. ,0.2,5) 
CALL GRID (3,-20. ,0. ,0.2,5) 
I A=-l 

DO 701 1=1,3 
YS=I A- ( 3-1 ) xO . 1 


CALL CHAR (-50 
WRITE (20,700) 
FORMAT (12) 


. , YS, 0.1,0. 1,0. ) 
IA 


IA=IA+1 


CALL PLOT (1,0. ,0.) 
DO 710 1=1 , 1025 
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X=I 

CALL PLOT ( 2»X* CTC 1 3 3 
710 CONTINUE 

CALL GRID (1, 1025. ,0.,!., 13 
CALL GRID C3,1025.,0.,l.,13 
CALL CHAR ( 1050 . , 0 . , 0 . 25 , 0 . 25 , 0 . 3 
WRITE (20,7203 NCL 

720 FORMAT ('(',12,* 3*3 

CALL GRID (2,1000. ,0., 100. ,103 
IF (NCL . LT. 63 GO TO 730 
CALL CHAR ( 1050 . , -1 . 75 , 0 . 1 , 0 . 1 , 0 . 3 
WRITE (20,7213 

721 FORMAT ( ’MS’ 3 

CALL GRID (2, 1000. ,-1.5, 100. ,103 
TIME =0 
XI- -50 
DO 725 J = 1 » 6 

CALL CHAR (XI , -1 . 75 , 0 . 1 , 0 . 1 , 0 . 3 
WRITE (20,7233 TIME 
723 FORMAT (F5.03 

TIME = TIMF +100 
725 X1=X1+200 

CALL PLOT (1,1500. ,0.3 
GO TO 750 

730 CALL PLOT (1,0. ,0.3 
750 CONTINUE 
RETURN 
END 
C 
C 

SUBROUTINE PDIPO(L, FO ,NEQ , DE , NP3 
C 
C 

C THIS SUBROUTINE PLOTS THE GEOMAGNETIC FIELD LINE THROUGH WHICH THE 
C VLF PULSES PROPAGATE 
C 
C 

REAL L»FO,NEQ,XC( 4 9 3 , YC( A9 3 , XD( 21 3 , YD( 213 , ANC 1 1 3 
INTEGER DE , NP 
PAI=3 . 16159 
C CALCULATE A CIRCLE 
DO 600 1=1,26 
XC (1 3 = C0S(PAIK(I-1 3/263 

600 YC ( 1 3 = SIN(PAlX(I-l)/26) 

DO 601 1=25,68 

11=1-26 

XC ( 1 3 =-XC ( 111 

601 YC ( 1 1 = -YC ( 111 
XC ( 69 3 = 1 

YC ( 6 9 1 = 0 

C CALCULATE A DIPOLE FIELD LINE 
RT=SQRT(6600»L/6500-ll 
FAI1=ATAN(RT3 
DFP=FAII/10 
PA =-FAI 1 
DO 610 1=1,11 
SI=S1N( PA3 
Cl =C0S ( PA 3 
AN(I3=PA«180/PAI 
R=L»CI*X2 
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PA=PA+DFP 
XD(I)=R*CI 
61 D YD C I ) =R*SI 

DO 611 1=12,21 
16=22-1 
XD ( I ) =XD ( 16 ) 

611 YD( I ) =-YD ( 16 ) 

CALL RAML 

CALL SCALE Cl.,1.,-1., 2.) 

CALL PLOT (-2,1. ,0.) 

DO 620 1=1,69 

620 CALL PLOT ( 0 , XC ( I) , YC C 1 3) 

CALL PLOT ( 1 , XDC 1 ) , YDC 1 ) ) 

DO 630 1=2,21 

630 CALL PLOT ( 2 , XD ( I ) , YD ( I ) ) 

DO 632 1=1,11,2 

CALL PLOT (1,XD(I),YD(I)) 

CALL POINT! 2) 

1S=( 1+1 )/2 

X1=XD(I 1+15*0.06-0.2 
Yl=YD(I)-0. 5+15*0. 05 
CALL CHAR ( XI , Y1 , 0 . 1 5 , 0 . 15 , 0 , ) 

WRITE (20,631) 15 

631 FORMAT ( ' ( ' ,11, ' ) ' ) 

632 CONTINUE 

CALL GRID(G ,0. ,0 . ,1. , 5) 

CALL CHAR(-1. »-3. ,0. 15,0.15,0. ) 

WRITE (20,633) L.DE.NEQ 

633 FORMAT! ' L= » , F5 . 2 , 3X, • DE » , II , • MODEL * , 2X , ’ NEQ= ',F7.2, 

1 ' PER C.C. ' ) 

YD1 =-3 . 5 

DO 660 1=1,6 
16=1*2-1 
AG=AN( 16 ) 

CALL CHARC-1 . ,YD1, 0.15, 0,15,0.) 

WRITE (20,635) I, AG 

635 FORMAT CLAT. OF LOCATION ’,11,’ IS ' , F6 . 2 , ' DEG*) 

660 YD1 =YD!-0 . 5 

CALL CHAR (-1. ,-6.5,0. 15,0.15,0. ) 

WRITE (20,661) FO.NP 

661 FORMAT ('CARRIER FREQ, = ',F6.1,' H2',2X,' PULSE LENGTH= 

2 1 3 , ’ MS ' ) 

RETURN 

END 
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Appendix B 


VLF WAVE-PARTICLE INTERACTION IN THE MAGNETOSPHERE 

In this appendix, we shall discuss VLF wave particle interaction 
(WPI) processes in the magnetosphere. It is not our intention to give 
a complete theory for VLF wave amplification and triggering but, rather, 
to give the background and rationale for the VLF wave injection experi- 
ments related to this repox’t. 

The WM wave is assumed to propagate along a static magnetic field 

line. The electrons travel in the opposite direction. The motion of an 

electron in a WM wave embedded in a uniform magnetoplasma is studied 

first. We shall use a simple model to illustrate phase bunching and 

phase trapping mechanisms, as well as the relationship between bunching 

time and coherence bandwidth. Phase bunched currents due to q(V XB ) 

X w 

forces will be illustrated schematically. Then, electron motion in a WM 

wave in the magnetosphere is discussed. We then discuss debunching mech- 
anisms. 

According to the Lorentz force law, the electron motion in a TO wave 
embedded in a homogeneous magnetoplasma is 

mV>= ?= -eCE^ + V X B y ) (B.l) 


This can be written as [e.g., Bell, 1965; Dysthe, 1971; Matsumoto, 19723 

eB 

\ - L ir siI1 * ■ 

\ ■ '( v « + f)-ir sin * te.2 b ) 
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where 


/ \ eB 

* - fl o - k ( v „ - e) - TT 


to 

eB V + r 

W I! k 


cos i}r 


(B .2c) 


ft = eB /m: electron angular gyrofrequency. 

o o 

to — WM wave angular frequency 
1c = wave number’ 

\]r = angle between -B . and 


We have used the relationship 


E 


w 


B 


w 


to 

k 


(B .3) 


The geometry of B q , V^, V , k, B^, and E w is shown in Pig. B.l. 

From Eq. (B.2), we can obtain two constants of motion [Bell, 1965; 


Dysthe, 1971; Matsumoto, 1972] 



where 


(B .4a) 


(B .4b) 


ft .» 
w 


eB 
w 

m 


ft - to 

V s- — — is called the cyclotron resonance velocity 

o k 


Equation (B. 4a) implies that electron energy is conserved in the 
frame moving with the wave phase velocity to/lc. It is important to em- 
phasize that wave and electron are traveling in the opposite directions. 
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/ 

/ 
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Fig. B.l. CONFIGURATION OF E AND B FIELDS OF A 
WHISTLER MODE WAVE AND VELOCITY COMPONENTS V 
AND V x OF AN ELECTRON IN THE LAB FRAME. The 
wave and the electron ti'avel in the a 2 and -a 2 
directions, respectively. \[r is the angle between 
~Bty and V^ . B 0 , a constant magnetic field, 
aligns in the £ direction. 
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In that frame, there exists only static magnetic fields. The forces ap~ 
plied to electrons are always in directions perpendicular to the electron 
motions . There is no energy exchange between the wave and the electrons . 

The second constant of motion X has been called the "inner Hamil- 
tonian" [Liemohn and Duane, 1976]. 

In the 3-D phase space where V ^ , and i|r are chosen as the 

2 

three orthogonal coordinates, the two constants, U and X, represent 

2 

two surfaces . A constant U describes a cylindrical surface of radius 
U, where U > 0, centered at V = co/k and V = 0. The axis is par- 
allel to the i|r-direction . 

A constant X represents a rather complex surface in the phase 
plane. The cross section of the surface in a plane with constant ijr is 
a parabola with the axis parallel to the axis. The vertex is at 

CV , -(I«X/2ft cos f), i]f] , and the length of the latus rectum, the chord 
drawn through the focus and perpendicular to the axis of the parabola, 
equals 2fi^ [cos \|rj. The parabola is open towards the positive V di- 
rection when -it/2 < i]r < n/2 and towards the negative V direction 
otherwise. We are interested in the portions of the surface on which 
V > 0. Figure B.2 illustrates the loci of the vertices of the parabolas 
in the phase plane for three cases : X > 0 , X = 0 , and X < 0 . The cor- 
responding surfaces are shorn in Fig. B.3. 

Given the values of V , , and i|r of an electron at an instant, 

2 

U and X are defined thereafter. The trajectory of the electron in 

phase space is the intersection of the two surfaces defined by the values 
2 

of U and X. Figure B.4 shows several electron trajectories in which 

we use the same V , V but various W for the electrons at an instant 
x II 

2 

to calculate the values of U and X of the electrons . It turns out 
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Fig. B.4. TRAPPING CONTOURS IN THE 3-D PHASE SPACE. 
0 e is the small angle such that electrons with 

v n = v o at If I > it - 9 £ can not be trapped by 
the wave . 
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that the value of U is the same but that X has various values . 
There are two kinds of trajectories. Those with closed contours corre- 
spond to the trapped electrons . The others with open contours are the 
trajectories of untrapped electrons. There are -several important fea- 
tures worth mentioning. 

(1) All the trajectories are on the same cylindrical surface 

2 

defined by the first constant, U . 


(2) Not all the electrons whose trajectories pass through 

V = V q can be trapped by the wave. Those which pass 

through V ([ « V q at [i|r| > jt-0^ are not trapped. An 

approximate expression for 9 has been worked out by 

Gendrin [1074] . We have independently worked out this 

angle numerically for a model in which f n = 100 kHz , 

2 y 9 

f fl * 9.8 kHz, f - 4.5 kHz, and U = 7 xl0 14: Cm/sec) . 
The result is shown in Pig. B.5. For B <10 my, this 
angle is very small (<1°). However, for waves with 
B ~ 100 my, 9 can be as large as 5°. 


There are two important points in phase space: the so-called resonance 

point at 


fl / fl /kV > 

ilr = 0 and V = V ^ V ( W - -A — 

ii r - o k \1 + fl /kV 
\ w J.; 


(B.5a) 


and the anti-resonance point at 


fl / fl /kV N 

t = ±.Jt and V = V ^ V + L W n .TV 

it ar o k \1 - fl /kV 

\ w 


(B.Sb) 


They are the solutions of Ecl, (B.2) when V. = = Tjr = 0. Furthermore, 

the resonance point is a stable point and the anti-resonance point is 
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(degree) 



B w (mr) 


PITCH ASS ™ G EIJ3CTR0N 

greater than 4° when the wave is strongei^than ^OcTmy. 



generally an unstable point. The projection of these points, as well as 

that of electron trajectories on a V ( plane, are shorn in Pig, B.6. 

It is interesting to point out that an electron which passes through 

V n “ V ar’ re £ ardl ess of the value of i{r, is always trapped by the wave. 

The deviation of from V is proportional to B for small B . 

w w 


At two extreme conditions, the two resonance velocities V and V 


are identical. When Fl / k » V 

w x' 


ar 


V = V = 
r ar 


w _ 

k ~ “ V P h 


CB.6) 


By recollecting that V ( and Vph in our convention are in the opposite 
diiections, we realize that Eq. (B.6) is a longitudinal resonance condi- 
tion similar to (but not) the Landau resonance. Whe n V » Q /k 

i w ’ 


V = V — V 
r ar o 


(B.7) 


Under this condition, V q is the "real" cyclotron resonance velocity. 

The trapping width in V (j is not symmetrical about V = V The 

lower boundary extends further in V (| than the upper boundary. But the 

difference is very small. Figure B.7 shows the difference AV - AV 

tL tU 

ana the total width AV tL + ^ for a Model in which f p = 100 kHz , 

f H = 9 ' 3 ltHz ’ f = 4 - 5 kHZ > electron pitch angle a » 35”, and u 2 = 
14 2 P 

\ x0 (m/sec) . When B w < 10 my, the ratio of the difference to the 

total range is less than 1 percent. 

A way to obtain an approximate formula for the trapping width is to 
assume a constant V ± instead of using the first constant of motion U 2 . 
From Eq. (B.4b), we can find the value of X for the electrons whose tra- 
jectories are on the boundary of the trapping width by setting V = V 

II o 

at ijr = +*, i.e., 
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Pi S; *il: bTo ™T(v 0K -°* f , t ™ d T pibg C0!TODHS shom 

pass " 04 AV “ ^ 

I"v = . T ° ”* ♦ ■ “ 

II v o 
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Fig. B.7. THE TOTAL TRAPPING WIDTH AV t y + AV tL AND THE 
VALUE OF AV. bL - AV^ VS WAVE INTENSITY. 
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AV it _ + AV tu (10 m/sec) 



(B.S) 


A \ 

X. = -A 

o wy k j 


The values of V at * = 0 under the condition x-^ are 

/I) V \ 1/2 ' 

V = V ± 2[- W x . 1 

Hot k J 


(B.9) 


The trapping width in v is therefore 


/n v \ 1//2 

iv t - «(T i ) 


(B .10) 


Usihg the same mode! as the one In Fig. b.T, we plot hot], A V) . of Eq. 

(B. 10 ) and the 'Wet" trapping width AT. L + AV tff on Mg. 3.8. It !s 

noticed that the discrepancy Is very small even for & w - 150 my. In 

general, Eq. (B.lo) is a good approximation of the trapping width fox- 

electrons with pitch angle greater than 10 degrees and wave intensities 
less than 150 1117 . 

So far, we have discussed only electron motion in a TO wave in a 
Phase space in which the temporal variations are not explicitly shown. 

It would be beneficial to illustrate the electron motions in a TO wave 
in the time domain. 

We shall start with a brief discussion of the importance of the terms 

on the right-hand side (bus) of Bq. (b. 2). T1 ,e first terms in Eqs . (3.2a) 

and (3.2b) are due to V ± », and V. X5< forces, respectively, me 

second term in Eq. (a.2b), (- uA ) (eB/m) sin *, is due to the E force. 

w 

In Eq. (B. 2c), the first term J 1 q and the second term k(V +wA) rep- 
resent the V ± X B q (gyration) force and the doppler shift to the wave, 
respectiveiy. The third term comes from \X\ and E^. forces . 
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AV x 10 m/sec 



Fig-. B.6. THE CALCULATED TRAPPING WIDTW Av ov , 
AM THE VALUE OP AV + Av “ “ " (B 


VS WAVE INTENSITY. 



Equation (B.2c) can be written as 


t 


k(V - V ) 
o || 


eB V 4. - 
w li k 
U V — GOS * 


(B.ll) 


liie last teim is, in fact, the very term that causes V and V to 
. r ar 

eviate from The effect of this tern on the trapping region becomes 

pronounced near * = 0 and ±*. We have show earlier that this 
effect is relatively small for B, < 10 m 7 and electron pitch angle >10°. 
This term will become important when electron pitch angle is small [Bye - 
the, 1971; Inan, 1977] even in a weak field. 

By neglecting the last term in Eq. (B.ll), W e have 


t - MV - V ) 
o II 


We can lump Eqs . (B.2a) and (B.12) together into 


(B .12) 


/kV eB \ 

^ + (' — ~ — sin ijr = 0 


in 


(B.13) 


This is a pendulum equation with B w analogous to the gravitation force 


g and with m/Cekv^) analogous to the length of the pendulum f. A 
change in V ± implies a variation in f. per the purpose of illustra- 
tion of the XB w phase bunching mechanisms, let us assume a con- 
stant V A - This analogy is illustrated in Pig. B. 3 . The period of this 

oscillation for a small angle approximation is T , where 

t 


T = 
t 


2jt 


eB 


1/2 


(B . 14a ) 


kV. 


w 


l m 
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Fig. B.9. AN ANALOG OF AN ELECTRON MOTION IN A WM WAVE TO A PENDULUM 
MOTION IN THE GRAVITATIONAL FIELD. The time domain motions are shown 
in the upper panel, and the corresponding: phase domain motions are in 
the lower panels. 
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T is the trapping period for small phase angle electrons . 
t 

As the pendulum swings from Locations 1, 2, 3, to 4, it executes a 
trajectory in the phase plane (ijr,\|r) as indicated by 1, 2, 3, and 4. 
Similarly, as an electron in a mi wave with a phase angle t changes 
from positions 1, 2, 3, to 4 in the time domain, its position in the 
phase plane (V ( ,\jr) is moving from 1, 2, 3, to 4, as indicated in Fig. 
B.9. The condition V = V Q implies \jr = 0 and occurs only at the two 
locations where |i|rj = maximum. In this example, V ( = V q occurs only 
at Locations 1 and 3 . 

Without any other external force, the pendulum will exercise a pe- 
riodic (but not sinusoidal) oscillation in time if the condition = 0 
occurs within -130° < %jr < 180°. Otherwise, it will rotate. Similarly, 
the electron is oscillating or rotating if the condition V =V q can or 
cannot be met within -180° < \Jr < 180°. The electron oscillating in time 
has a closed contour in the phase plane and hence is trapped by the WM 
wave. The one rotating in time has an open contour and hence is untrapped 
by the wave . 

For those trapped electrons with large phase angle i|/^ at V = V q , 

the trapping period depends on the value of i|r and is [e.g., Matsumoto, 

R 

1972] 


where 



(B.14b) 




K(x) 1. the complete elliptic integral of the first hind. K(sl n t /2) , 
as a function of * R , is plotted in Fig. B.10. I*,r electrons with\ R = 
-if at V t - V Q , the trapping period is infinite. 

^ shall recall that the period of a pendulum is almost constant at 
the value 2n ■JITi for reasonably small angles * Csay,15»). aiB is 

the hey to understanding the phase bunching mechanisms. Suppose there 

ate ^ny electrons initially at V, . V o hut at various phase angles. 

Their V/e are the same, initially, there is no transverse current 

as she™ in Fig. B.ll. Corresponding to the analogy of pendulums, the 

pendulums are identical and are at various locations but at rest (^,= 0) 

initially. As time goes on, all the pendulums converge towards * = o. 

At one quarter of the period, n/S Vl/g, the pendulums initially a- small 

* with each other at * = 0. m the wave-electron 

system, this phase angle convergence mechanism is .ailed phase bunching. 

It is this process that creates a transverse current in the direction of 

B w , as shorn in Fig. B.ll. me bunching time T b is defined as a quar- 
ter of T. . 


1/2 




b 2 1 eB kV 
\ W 1/ 


(B .15) 


It is interesting to show that nc + 1 , 

o 15 the corresponding frequency 

range of the trapping width, which we shall call +i, , 

, on we snail call the coherence bandwidth 


Af ^ .Ji, Aar _ ^ - .1/2 — [ 

CB = 2* AV t ~ * (kV . a J ~ T, 


1 w 


(B .15) 


It is important to point out that these electrons at the same loca- 
tion in space initially will not arrive at another location later simul- 
taneously because their V, 's are not identical (due to different v^XB 


Fig. B.10. 

vs V 


PLOT OF THE VALUE OF 


K( S in ir / 2) 

U 


IN EQ. (B .14b) 



Fig. B.ll. 


ILLUSTRATION OF PHASE 


BUNCHING MECHANISM. 


See text. 



forces) [Helliwell and Crystal, 1973). There is a spatial spreading in 
the phase-bunched current. In other words, a single sheet of electrons 
at the same location with the same and V but various i|r initially 

can be phase— bunched by the V X forces to create a transverse cuz’- 
Z’ent which is in alignment with B^. The transvei’se cuzu’ent is not lo- 
calized but spread in space because the V 's of the electrons have been 
altered differently by the bunching forces. 

For a general case, the tz-ansverse current produced by the phase- 
bunching pi’ocesses can be decomposed into two components : parallel 

to the B w and J 2 perpendicular to it, as indicated in Fig. B.12. 


T U,t) = f f f qf(x,V ,V ,t,t) V cos t dijr dV dV (B.l7a) 
Jy Jy 11 ■*- x X ll 

If X 


J 2 (x,t) = f f f qf(x,V ,V ( i|r,t) V 

Jy Jy J-.fi ll -*• 


x sin i|r dzjr dV^ dV (B.17b) 


Even though there is no general agreement on how large the ranges of V 

il 

and must be over which the integrations shall be carried out [Sudan 

and Ott, 1971; Dysthe, 1971; Nunn, 1974; Crystal, 1974], the latest com- 
puter simulation results [Cz*ystal, 1977] seem to confirm that integrating 
only over the trapping width can give a good indication of the magnitude 
of these simulated currents.- Furthermore , it is shown in Nunn [1971] 
that, in an inhomogeneous medium, "cyclotz’on resonant electrons may become 
stably trapped, undergoing a steady change in energy and momentum. Where 
wave-particle-intez’action time exceeds one trapping period, bunches of 
stably ti'apped particles completely dominate interaction with the wave.” 
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We believe that both trapped and nearly-trapped electrons are im- 
portant to the wave growth. However, since trapping is such a simple 
idea that can give us good indications of the stimulated currents and, 
since we are interested in a qualitative answer at the moment, we shall 
consider the trapped electrons only from now on. 

As to how these currents contribute to wave amplification, there 
are different opinions, too. In the model of Helliwell and Crystal 
[1973] , the current is treated as a source current which radiates 

a new field. The new field combines with the old wave, resulting in net 
wave growth. Frequency variations come from the inliomogeneity in the 
geomagnetic field lines. In Dysthe [1971] and Nunn [1974], a gradient 
in electron V distributions is required to produce J . Since <T 

II 2 2 

is aligned with E^, there is energy exchange (J 2 • E ) between the 
wave and the electrons . serves as the source of frequency variations 

in both homogeneous and inhomogeneous cases . 

Hie inhomogeneity of the geomagnetic field in the magnetosphere is 
important in WPI processes . For instance , the inhomogeneity acts on the 
electrons like an external torque, pushing the stable point in -t|r from 
- 0° to other values, deforming electron trajectories in phase spaces, 
and reducing the trapping width. Far away from the equator, the inhomo- 
geneity is so large that trapping becomes impossible [Dysthe, 1971; Nunn, 
1974] for CW waves with reasonable amplitudes. More quantitative argu- 
ments on when tripping is possible will be discussed later in this sec- 
tion. 

The inhomogeneity forces in the equations of motion are modeled 
[e.g., Dysthe, 1971] as 
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mV 

F = — i 
II 2B 5z 


(B.lSa) 


„ mV V ^ 

p - 1 H oB 

II 2B Sz (B.18b) 

_ They "* derlVad from fi-t adiabatic constant and the Gauss law 
V -S = 0. it has bean assumed that the spatial variations in B are 
very small over a distance of the order of the electron gyro-radius. The 
forces have been averaged over a gyro-radius. This is a widely-used 
model [Dysthe, 1971 j Nunn, 1971; man, 19773. 

Electron motion in a W wave in the magnetosphere, therefore, can 
lie modeled as [e.g., Dysthe, 1970] 


. eB V N 

V. = V — - sin \(f + -A 

Y 2b 5z 


± m 


CB.19a) 


where 


V = 


w 


i eB 


l w 

sin ijr 


+ r) 

eB 

w 


m 


W s n 
X II OB 

2B 3z 


+ k 


COS ijf 


= eB/m 

R E ~ earth r adius 
I* = Mcllwain's L parameter 


(B .19b) 


(B .19c) 


We have assumed that V - 0 is at the eguator and that the wave and 
electrons travel in « and -a directions, respectively. The curvature 
of the geomagnetic fields is not included in the formulation. The geomag- 
netic field near the equator can be assumed as 
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2 K 


* 

a 

z 


(B .20) 


Numerical solutions of these 
1977J for various values of b . 


equations have been calculated Cinan, 


Equation (B.19c) is the same as En m . 

S Eq - (E - 2o) - By the same argument 

bef oi'e , f or electrons with a pitch angle „ _ 

s gieatei than, say, 10 de- 

"* 0 P0BdUlU " 1. Obtained as iollows: 

f « n - k hr + iiA 

V' 1 V (B .21) 

By differentiating Eq. (B.21), We have 


iff = Q - kV - kV 


» ll 


We have assumed a constant frequency signal so that = Q 
«.e UK dispersion relation, it can be sho lm that 


(B .22) 
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ubstituting CB ,19a) and (B. 23) into (B.22), we have 
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Let us define the RHS of (B.24) as r/m, i.e., 


r/m £ - 


W* 

X 

2(1 


+ V 


kV \ 

» 

2(G - u)i 




(B .27) 


The inhomogeneity works as an external torque t which is a function of 
position. When |t| < eB^cV^ , we can define an angle ijr Q so that 

sin = T/eB lcV (B.28) 

T o w JL 


Equation (2.22) can then be written as 


.. kV eB , „ v 

i|r + i — - [sin - sin ijr 3 « 0 (El. 29) 

At ijr = ijr , ijr = 0. The inhomogeneity force has shifted the stable point 

of the pendulum from * = 0 to . Electrons are still trapped by the 

wave but they oscillate about i|r = i|f o instead of about i|f = 0 when T 

varies sufficiently slowly for a distance of the order of V |( • T t> where 

T is the trapping period. The trapping width is reduced, too [Dysthe, 
t 

1971; Nunn, 1974]. On the other hand, when |t | > eB w kv x » trapping is 
impossible . 

At the end of the triggering wave, assuming that the stimulated cur- 
rents are not strong enough to maintain a self-sustained system, the or- 
ganized electrons will undergo a phase-mixing process which reduces 
transverse currents. In a homogeneous model, this process occurs because 
of a finite V spreading in the organized electrons. Let us assume that 
only trapped electrons are important in the bunching processes . At the 
end of a wave train where t = 0, B w disappears. From Eq. (B.13) for 

t > 0, 
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(B .30) 


i|r = 0 


Integrating twice, we have 


t ~ At + B 


(B .31) 


where A, B are constants to be determined by the initial values of ijr 
and Tjr at t = 0. 

From Eq. (B.12), 

i=k(V o -V (i ) (B .32) 

Thus , 

t = K(V q - V (| ) t + B (B .33) 


The trapping range in the phase plane is deformed, as shown in Fig. 
B.13. Note, the V ( ’s do not change but the ljr's do. The electrons 
eventually will be totally phase-mixed. The drift in ijr is ^ . linear 
function <V Q -V |( )t. The further away from v , the larger the drift 
in i|f. The largest difference in ty-drift is between the electrons with 
V„ = ± AV t /2 and can be written as 

^(t) = -2jrAE CB t (B .34) 

The debunching time is defined as the time required to randomize the or- 
ganized electron. It can be shown, when ^(t) = ±2jt, the transverse 
currents totally disappear in this case. Thu%, the debunching time T D 
can be obtained as 


T = — - 
D Af 


- m 


CB 


(B .35) 


In an inhomogeneous model, in addition to this finite spread of ?'s 

II 

that causes phase-mixing, there is another adiabatic force which helps 


169 



Fig. B.12. CONFIGURATION OF Jg OF 

EQS. (B .17a) AND (B.17b). 



Fig. B .13 . VARIATION OF A TRAPPED REGION AT THE END OF 
TRIGGERING WAVE IN A HOMOGENEOUS MAGNETOPLASMA. 
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the phase mixing. As can he seen from Eq. (B.24), various values of 
either or V will produce different external torque t* This 

causes differential drifts in ijr to produce phase-mixing. Thus, in 


an inhomogeneous media, the debunching time should be smaller than T . 

D 


Appendix C 


DERIVATION OP SIDEBAND REPRESENTATION OF FSK WAVES 


In this appendix, we shall derive the sideband representations of 

the FSK waves analytically. The sideband amplitudes, frequencies, and 

phases will be given quantitatively. 

An FSK wave (except the calibration wave) consists of a sequence of 

RF pulses , each lasting r msec , alternating between two carriers sep- 

S' 

arated by &£ , It can be viewed as two pulse trains: one at the car- 
rier f^ and the other at f 2 = f^ -Af. Each of the pulse trains is a 
square-wave modulated sine wave consisting of a sequence of RF pulses 
and gaps in between. All the pulses and the gaps have the same duration 
T g • The phase difference between two successive pulses £0^ equals 
2jtAf T g which creates either an in-phase or an anti-phase condition. 
Analytically, the FSK wave can be written as 

= °n( t 

+ f3M (t) sin (2jtf t + 0. + 0)1 (C.l) 

O y n c* 1 O J 


)[' 


M (t) sin (2nf.t + 0 ) 
o,x 1 o 


where M (t) are the square wave trains shown in Fig. C.l. The sub- 

O j3t 

scripts o and n denote the 0- and ^-conditions , respectively. M (t) 

O , 7t 

extends from t = -oo to +oo. c and CS are the wave intensities of f 

K 1 

and f p , respectively. The inherent limitations on the transmitting 

system at Siple cause [3 to be different from unity and 3 < 1. The 

1 aviation usually is small when Af « f . f"](t) is a square pulse with 

a unit height and one-s duration centered at t = 0 . 0 q is the initial 

phase angle of the carrier at f.. and is constant. 0. equals 2itAfT • 

X s 


M 
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Flg * C,:L ' PLOT OF THE FUNCTTnwc; M r^s 

EQ. (C.2b) VS TIME. T il 2? V*) OF 

t s xs the small pulse length, 


Without losing generality, we can assume 0 = o and C = 1. 
Eq. (C.l) can be written as 


Then , 


W 


= n(* w 0 (jt (t) 


(C.2a) 


where 


V* (t) - ■»“ **V * (*,„<* - t.) sin ( 2 « 2 t te . 2b) 

The Fourier transform of W Q ^(t) is defined as 


= J [ W o,, (t >] = jf W D> /t) „ 


(C.3a) 
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■ r v„<« •* 

“'—CO ? 


i2jfft 
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A o,« (f) can 8e written as Ce .g . , Bracewell, 1965] 


A (f) = m 
o,Jt o 


,* a> * si [ B « - v - bo; + fi >] 


“i2jtfr 


+ /3e 


m o,« (f) * £ 2i~[ 8(f "V “ 5Cf + V] <C.4) 


where 


m ^ Cf) 
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e dt 


SCO = Dirac delta function 


* = convolution 
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are 


The amplitudes and phases of m (f) and ra e “ l2ltf T s 

o,jr 

shown in Figs. C.2 and C.3. They are the base-band spectrum of W (t) 

o,jt 

consisting of discrete spikes separated at multiples of (2 t )“ 1 . For 

s 

example, the spacings are at multiples of 50 Hz when T a 10 msec. 

s 

The delta functions shift the whole base-band spectrum to the ap- 
propriate frequency bands centered at ±f and ±f . i n general f 
and f 2 are in kHz ranges; the spikes shifted from the base-band toward 
the positive, and those towards the negative frequency ranges do not in- 
terfere with each other. Furthermore, because the transmitted signals 
are real signals, the part of the spectrum on the positive and that on 
the negative frequencies are conjugate to each other. 

That is , 

A(f) = A + Cf) + A_(f) ( C .5a) 

and 

A_(f) = A*(-f) (C. 5b) 

where 


, % 1 -i2jtf<r 

■ m * 2l - V + 


h -± 0 . 


(C.5c) 


Figure G.4 shows the amplitudes and phases of A + (f) for the fol- 
lowing two cases: (a) <r g . = 10 msec and Af =100 Hz, and (b) t =20 

ms and Af = 100 Hz. Both are 0-waves. Figure C.5 shows the following 
two; (a) T g = 10 ms and Af = 150 Hz , and(b) r s = 20 ms and AO = 
125 Hz . Both are jt-waves . p has been assumed to be unity for simplic- 
ity. The frequency spacing between the spikes for the signals with t = 

s 

10 and 20 ms are multiples of 50 and 25 Hz, respectively. The spectrum 
A + (f) is symmetrical about the frequency at (jL +f„)/2. 

JL *d 
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It is noticed that there are no spikes at ^ and f g , the carrier 
frequencies, when the signals are rt-waves . This is understandable. The 
odd multiples of *-radian phase shifts between successive pulses in a 
carrier have cancelled the signal intensity completely in the carrier 
frequencies, it is a kind of carrier suppression transmission. There- 
fore, there is no energy at the carrier frequencies. 

A o,* (f) consists °f discrete spikes in frequency. From Eq. (C.3b), 


we have 


* T 


Pj sin (2«f t + 0 ) 


(C.S) 


Each spike corresponds to a sinusoidal wave in time domain. The FSK wave 


W o 7 r (t) can be represented as 


o,jr (t5 “ n ~ | j J Pj sin (2jtf jt + 0 ^ 


(C.7) 


It consists of several sinusoidal waves at various frequencies. Each 
of the sinusoidal waves last for one second, exactly the same duration as 

the FSKed wave. These sinusoidal waves , or spikes in the spectrum, shall 
be called sidebands. 
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Appendix D 


MISCELLANEOUS EXAMPLES 

The purpose of this section is to show a collection of interesting 
VLP wave phenomena related to "coherence bandwidth." We shall present 
some of the experimental results from two of the VLF wave injection ex- 
periments . 

1. COBA Program 

There are two kinds of signals in the coherence bandwidth (COBA) 

program. The first kind involves stepping up or down in frequency in a 

one-s pulse at the Tth msec. The frequency offset Af ranges from -25 

to 25 Hz . The first r msec of the pulse organizes electrons in a range 

AV centered on the corresponding resonance velocity. When the fre- 
li 

quency offset Af is greater than the coherence bandwidth, most of the 

electrons organized by the pre-offset section of pulse become "free." 

They may radiate a field in a natural mode. Thus, triggered emissions 

may be developed at the offset. On the other hand, when AC is less 

than the coherence bandwidth, most of the electrons organized by the 

pre-offset section can interact with the post-offset section. We may 

only expect to see a small amplitude perturbation at the transition. 

The second kind is the so-called "rake-spectrum" signal. A pulse 

at f is transmitted for 400 msec, followed by an FSK wave with t = 
o s 

10 ms alternating between f ± 75 Hz for 1.1 sec. The rake-spectrum 
signal is illustrated schematically in Fig. D.l for both the high time 
resolution and the high frequency resolution displays. The intensity 
of the sideband at f is about one-third that of the pulse. The phase 




kHz 



-J 1 L. 


FI ~ at a high ™» 



is not continuous at f, during the trausitioh. Thera is a difference 
of 3g/4 radians between tie phase of the hey down pulse and that of 
the sideband at t g . It is important to point out that the transition 
in the rake-spectrum is not smooth, because an electron can not resolre 
the sidebands in the interaction region until it encounters more than 
four 10 msec pulses. Thus, the key down pulse may trigger emissions 
that may then be entrained by the sidebands. 

This program produces a variety of interesting phenomena, but it 
has not been possible to provide conclusive evidence regarding the size 
of the coherence bandwidth because of the lack of a large set of data 
samples. Two of the most interesting phenomena occur around the vernal 
equinox and will be discussed in this section. 


Multiple triggering of rising emissions as shown in Fig. D.2 are 
among the typical features of wave activity during the equinox periods. 
The triggering signal is a one-s long pulse. It triggers a rising emis- 
sion at about the 420th msec. The emission starts to develop around the 
300th msec. We shall call this kind of emission "pre-termination emis- 
sions." The pulse triggers another rising emission at the end. 

Figure D.3 shows four one-second waves in the COBA program. The 
frequency offset in each wave occurs at the 400th msec. The offset fre- 
quency £f is +5, -5, +25, and -25 Hz, as indicated. The ramps indicate 
that there are two paths and that the one with the longer delay is the 
stronger. The ramps also serve as a barrier to prevent the emission 
triggered by one wave from contaminating the following wave. The times 
of offset are indicated by vertical bars on the top of the panel; the 
short bar for the first path and the long one for the second (a stronger) 
path. It is observed that the negative offset in frequency has inhibited 
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Fig. D .2 . AN EXAMPLE OF MULTIPLE TRIGGERING; A 
TYPICAL FEATURE OF TRIGGERING AT EQUINOX. 


RO 16 MARCH 77 




A f 5 Hz -5 Hz 


+ 25 Hz -25 Hz 


Fig. D.3. DEMONSTRATION OF A PART OF THE COBA PROGRAM TRANSMITTED 
AROUND EQUINOX. The pre-termination emissions are observable when 
the triggering waves have a positive frequency offset and disappear 
when the waves have a negative frequency offset. 
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or delayed the pre-termination emissions, while the positive offset in 
frequency does not affect the development of this emission. This indi- 
cates that a 25 Hz offset in frequency of a signal is not large enough 
to allow the two sections before and following the transition to become 
independent. Thus, 25 Hz is less than the coherence bandwidth of the 
signal in the pre-offset section in the interaction region. 

The mechanisms that change the form of the pre-termination emis- 
sions due to negative offset in frequency have not been identified yet. 
However, it is believed that the timing of the offset is an important 
factor, in all the cases shown in Pig. D.3, the offset time is at the 
400th msec, the time when the rising emission would have been develop- 
ing, had there been no offset. A positive frequency offset appears to 
favor the development of a rising emission. A negative frequency off- 
set seems to inhibit the development of rising emission. 

Figure D.4 shows an example of "rake-spectrums There were two 
paths, but the second path is the more pronounced. The 400-msec long 
pulse has been amplified. It may have triggered a BLI type of emission 
rHelliwell, 197S] that appears to feed energy to the 2.45, 2.55, and 2.5 
kHz sidebands. The sideband amplitudes appear to decrease and finally 
reach the steady state levels. The total amplitude has dropped by about 
10 dB from the peak value to the steady-state value, as indicated by the 
amplitude chart. The transient period is about 250 msec. 

One more example is shown in Pig. D. 5. There is only one observ- 
able path in this period. The emission triggered by the initial constant 
frequency pulse is entrained by the 2.55 kHz sideband and appears to 
fluctuate in between the two nearby sidebands at 2.5 and 2.6 kHz. The 

intensity of the emission forms a plateau for about 300 msec and is then 
quenched to a steady state level. 
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Fig. D.5. ANOTHER EXAMPLE OF A " RAKE -SPECTRUM" WAVE. The emission 
triggered by the 400 ms pulse is entrained by the 2.55 kHz side- 
band and wiggles in between two nearby sidebands at 2.5 and 2.6 
kHz. The intensity is almost constant for 200 msec and then drops 
10 dB to the "normal" steady state level . 
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It is postulated that this Phenomenon is a suppression effect. *. 
S 0 Hz spacing between tte has imposea Umitations m 

the steady state levels of the amplified sidehands . When the amplitude 
ot a sidehand is heyond its steady state level, it interacts with other 
sidebands at nearby frequencies. fbr reasons that we do not tally under- 
stand yet, the nearby sidebands have suppressed the sidehand with "ex- 
cess" energy. 

This postulate retires further investigation. It will be useful 

to transmit this program again near the equinox periods to collect more 
data samples. 

It is important to point out that the radiated energy from the trans- 
fer is distributed into the sidebands when the transmitter is operated 

ln a FSK mode. For example +h P 9 ^ it, 

■"Pie, the 2,55 kHz sideband has only about one- 

quarter of the total radiated energy, it shall not „ 

sy it: snail not be surprising to find 

coherence bandwidth of a single frequency signal is much greater 

than 50 Hz . In fact, the coherence bandwidth of the pulse near the tran- 
sition appears greater than 100 Hz 

Ite for the examples in both p igs . D.4 and 

D.5. 


2* CBT Program 

The coherence bandwidth trap (CBT) program is a complicated and long 
program. One cycle of the program lasts for four minutes and sixteen sec- 
onds. Its purpose is to study the effect on wave growth of shifting the 
frequency regularly * the ono-s triggering waves by a small amount; 
intend to present only a small portion of this program. In that portion, 
a signal is shifted between two frequencies f, and f, -« „ regular 
intervals of 10 msec. The phase is continuous during the transitions. 

190 


The frequency shifts are 1, 5, 10, 20, 30, 50, 100, and 150 Hz for vari 
ous signals. There is a one second "off” period between two successive 
signals in the frequency range (between f Q and f Q - 150 Hz). 

Figure D.6 illustrates the amplitudes of these waves of which f Q 
is at 3 kHz. It is obtained by band-passing the signals through a fil- 
ter at 3 kHz with a bandwidth of 170 Hz. The first two spikes correspond 
to the signals with Af = 0Hz, the second two to the signals with Af = 1 
Hz, the third two to the signals with Af = 5 Hz , and so on. It is ob- 
served that, as Af increases, the wave saturation levels drop. Figure 
D . 7 shows the plot of the saturation levels vs Af . Even a 1 Hz offset 
can cause a ~4 dB reduction in the saturation levels. This reduction 
becomes greater than 10 dB as Af increases beyond 50 Hz . 

The electrons cannot resolve the sidebands of these signals with Af 
less than 50 Hz in the interaction region because the modulation periods 
of these signals are greater than 100 ms . The electrons can only resolve 
segment-features of these signals. It can be shown that one hertz fre- 
quency difference can produce more than a jt radian phase difference be- 
tween two pure sine waves at ~3 kHz traveling from the ground through a 
duct at L = 4 to the equator. Thus, the propagation path between the 
transmitter site and the interaction region can introduce large amplitude 
distortions at the junctions of successive pulses, creating regular am- 
plitude depressions or enhancements before the signals arriving at the 
interaction region. The received signals show the effects of both the 
amplitude variation and frequency switching on the wave-particle interac- 
tion in the magnetosphere. 

Contrary to the previous cases, electrons can resolve the sidebands 
of the signals with Af = 50 , 100, and 150 Hz because of the relatively 
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Fig. D.7. THE PLOT OF THE AMPLITUDE VS M OF THE FSK WAVES 
ILLUSTRATED IN FIG. D.6. 
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small modulation periods on the signals (<40 ms). The propagation path 
introduces phase differences among the sidebands. Sidebands are constant 
frequency signals with a duration of one second (see Appendix C) . The 
absolute phase of a sideband is not important to its ability to interact 
with energetic electrons. Thus, the propagation distortion will have no 
effect on sideband growth. 

The sidebands are spaced at multiples of 50 Hz and extend in fre- 
quency more than the bandwidth of the filter (170 Hz), it is therefore 
essential to see the signal strength over a larger bandwidth. Figure D.8 
depicts the amplitudes of the signals over a 6S0 Hz bandwidth vs Af. 
Comparing this with Fig. D.7, it is noticed that the most significant 
change comes at the point of Af = 150 Hz. This is understandable be- 
cause the 170 Hz band -pass filter leaves out a significant amount of 
sidebands for the FSK wave with Af = 150 Hz . 

Figure D.9 shows the averaged wave intensities over four sets of 
data samples obtained within an hour. The bandwidth of the band-pass 
filter is set to 680 Hz. It appears that the wave intensity drops as 
Af increases for AC <100 Hz . 

For Af > 100 Hz, sidebands have been resolved by the electrons in 
the magnetosphere. Detailed studies reveal (not shown) that these side- 
bands are not amplified significantly and that there are no triggered 
emission. 

The 2.975 kHz sideband of the signal with Af = 50 Hz is about 1.4 
times stronger than that of the signal with Af = 150 Hz. The former 
can trigger emissions four out of five times, while the latter does not 
even show growth most of the time. This suggests that there is an am- 
plitude threshold level above which waves can be amplified and trigger 
emissions. 
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When Af < 50 Ha, there are effects in the data due to frequency 
offsets as well as amplitude variation introduced by propagation distor- 
tion. Frequency offset effects are small when A t is 1 He. Thus, a 
drop of a few dB caused by a 1 Ha offset most likely is due to the ampli- 
tude variations in the triggering signals introduced by the propagation 
distortions . 

Changing the frequency of a wave corresponds to altering the if of 
the electrons interacting with the waves Csee Appendix B) . Thus, some 
of the electrons trapped by a pulse at £ o may become untrapped by the 

following one at f, ±Af. But there are electrons which can effectively 

interact with signals at both frequencies. As Af increases, the number 
of these electrons which interact with signals at both frequencies is 
reduced. The ability of the wave to organise electrons is reduced and, 
hence, the growth is less. This effect should be small when the fre- 
quency change is small and shall become pronounced if the frequency shift 
increases. 

A way to eliminate the propagation distortion has been discussed in 
previous chapters • In order to study the effects on wave growth due 
to a small shift in frequency, it is necessary to reduce the propagation 
distortion by implementing a pre- and a post-processor on the transmit- 
ting and receiving sites, respectively. 
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